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We consider the local field dynamical temperature correlation function of the Quantum Non- 
linear Schrodinger equation with the finite coupling constant. This correlation function ad- 
mits a Fredholm determinant representation. The related operator-valued Riemann-Hilbert 
problem is used for analysing the leading term of the large time and long distance asymptotics 
of the correlation function. 

1 Introduction 

This work continues the study of the correlation functions of the Quantum Nonlinear Schrodinger 
equation (QNLS) out off free fermionic point which was originated in the papers 

We consider the temperature correlation function of the local fields of the QNLS equation, 



tr e T 

Here \l/(a;,t), \E'"l'(x,t), {x,t G R) are the canonical Bose fields obeying the standard equal-time 
commutation relations 



Abstract 




(1.1) 



[^ix,t),^^iy,t)] = 5ix-y) 



(1.2) 



and acting in the Fock space as 



vl/(x,t)|0) = 0, 



{Q\^\x,t) = 0. 



(1.3) 
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The evolution of the field \1/ with respect to time t is usual 

^(x, t) = e^^*^(x, 0)e-'"K (1.4) 

The Hamiltonian of the model H in ( |1.1|) and ( |1.4|) is equal to 

H = J rfx(9^^1'(x)9^^(x) + c^"f(x)^^(x)^(x)^(x) - /i^1'(x)^(x)) . (1.5) 

Here < c < oo is the coupling constant and h is the chemical potential. The parameter T in 
( |1.1| ) is a temperature. 

The quantum Nonlinear Schrodinger equation describes one-dimensional Bose gas with delta- 
function interactions. The basic thermodynamic equation of the model is the Yang- Yang equation 
1^ for the energy of an one-particle excitation e{X) in thermal equilibrium 

,(A) = A= - ft - ^ £ m (l + e-*) (1.6) 

The function e{\) behaves as — at A — > oo. It is positive for A G -R, if /?. < 0, and it has two 
real roots e{±q) = 0, if /i > 0. 

It is worth mentioning also the integral equation for the total spectral density of vacancies in 
the gas pt(A): 

/oo 2c 

where 

(1.8) 

is the Fermi weight. The value 'd{\)pt{\) defines the spectral density of particles in the gas. Due 
to the properties of e{\) the Fermi weight '(9(A) decays as exp{— A^/T} at A ^ oo. 

In this work we analyse a large time and long distance behavior of the correlation function 

The asymptotic evaluation of the correlation functions is one of the most challenging analytic 
problems in the theory of exactly solvable quantum field models. For the case of zero temperature 
the leading asymptotic term can be found via conformal field theory 0. The small temperature 
limit also can be considered in the framework of this approach, although, strictly speaking, increas- 
ing temperature destroys conformal properties of the model. In the present paper we develop the 
method which is based on the Fredholm determinant representations of the correlation functions 
and which allows to remove the small temperature restriction. 

A systematic exposition of the determinant representation method is given in |^. For the 
reader's convinience, we shall outline the principal features of the scheme together with a brief 
historical review concerning the asymptotic analysis of the correlation functions. 



1 + exp 



£(A) 

T 
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The determinant representation method is based on the remarkable fact that the correlation 
functions of the 1+1 exactly solvable quantum models can be represented as Fredholm determinants 
of the integral operators V acting in L^(S,(iA) and whose kernels, ^(Ai, A2), have the following 
special form : 



E;=ie,(A0/,(A2) 
Ai — A2 



5:e,(A)/,(A) = 0, 



:i.9) 



where functions ej(A), /,(A), and, in fact, the oriented contour of integration S, depend on the 
model under consideration. The first representation of this type was obtained in for equal-time 
correlation functions in one-dimensional impenetrable (c = 00) bosons. Later on, determinant 
formulae were derived for a majority of exactly solvable statistical mechanics and quantum field 
models (for the principal references we refer the reader to the book In particular, the Fred- 
holm determinant for the time- dependent correlations in one- dimensional impenetrable bosons 
was constructed in |10|. A generalization of the results of for non-free fermionic case was ob- 
tained in The determinant representation for the most general case of non - free fermionic 
time-dependent temperature correlation function ( |1.1| ) was found in 

The determinant formulae can be used to obtain nonlinear differential equations for quantum 
correlation functions. These nonlinear equations turn out to be classical integrable systems. Mor 
exactly, zero-temperature and equal-time two-point correlators are described by integrable ODEs 
of the Painleve type (see - flBI), while time-dependent and temperature or/and multi - point 
correlation functions appear to be the r-functions of integrable PDEs (see |TB[, [[I^] - ). It is 



interesting to notice (see e.g. |22|, and section || of this paper) that if the quantum system is a 
result of quantization of a classical integrable system, as is the case for the QNLS model, then the 
integrable PDEs describing the correlation functions belong to the underlying classical integrable 
hierarchy. 

The determinant representations for correlation functions can also be used to study their asymp- 
totics. In the zero-temperature case, a comprehensive asymptotic analysis of the correlation func- 
tions related to XXO and impenetrable Bose gas models has been carried out in [^, ||2^, |[T6[| , 
[ |17| , p5|] . For the two-dimensional Ising model the analogous results were obtained earlier in 



(see also |26| and 

A further development of the determinant representation method was achieved in the series 
PI 



of papers |2 



(see also g]). The approach of ^8 



is based on the use of the Riemann- 



Hilbert method of the theory of classical integrable systems for the asymptotic evaluation of the 
Fredholm determinants describing the correlation functions of the quantum integrable systems. 
The Riemann-Hilbert method allowed to extend the mentioned above zero-temperature results for 
the XXO and impenetrable Bose gas models to the general finite-temperature case (see ||28|, p9[| . 



and g). 

The Riemann-Hilbert asymptotic scheme is the principal analytic tool which is used in the 
present paper. In what follows, we describe its basic ideas in some details. 
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Riemann-Hilbert problem (RHP) appears in the theory of correlation functions due to a simple 
yet important fact that the resolvent kernel corresponding to kernel ( |1.9|) can be explicitly evaluated 
in terms of the solution of the matrix Riemann-Hilbert problem with the jump matrix G{X) given 
by the equation (cf. m, p|, [||): 

Gjk{\) = Sjk + 27iiej{\)fk{\). 

More exactly, let xW he a N x N matrix function which solves the following Riemann-Hilbert 
problem: 

1. xW ^ -^i ^ ^ oo, (normalization condition) (1-10) 

2. xW is analytic function of A if A ^ S, (l-H) 

3. x-(A) = x+(A)G'(A), AGS, (jump condition), (1.12) 
where x±W denote the (±) - boundary values of the function x(A) on S, i.e. 

X±(A) = hm y(A'), A'g(±)- side of E. 

A — >A 

Then, the resolvent kernel -R(Ai, A2) corresponding to the kernel (|1.9| ) {1 — R = (1 + V)~^) is given 
by the following explicit formulae, 

fl(Ai.A.)^ ^^-'f'Y'''^' . (1.13) 

Ai — A2 

N N 

E,{\) = E(x+(A)),.e,(A), F,(A) = ^ /,(A)(x;^(A))fc,. 

k=l k=l 

The dynamical parameters, i.e. distance x and time t, enter the jump matrix G{X) through 
the transformation, 

G(A) ^ e^(^'"'*)G(A)e~^(^'"'*), (1.14) 

where the rational in A and linear in x, t diagonal matrix function D{X; x, t) represents the disper- 
sion law of the underlying classical model (e.g. D{X]x,t) = diag(itA^ — ixX, —itX^ + ixX) for 
the NLS equation). It can be shown (see e.g. P]; see also [0) that the logariphmic derivatives 
of det(l + V) with respect to x and t (and, in fact, with respect to any other physical parame- 
ter) are explicitly expressible in terms of the resolvent R. Hence the asymptotic analysis of the 
original Fredholm determinat is reduced by (|1.13| ) to the asymptotic analysis of the oscillatory 
Riemann-Hilbert problem ( |1.10[ - |1.12| , |1.14| ). 

In the theory of classical integrable systems, the Riemann-Hilbert problems of the type ( |1.1U[ - 
|1.12| , |1.14|) represent solutions of the Cauchy problems for integrable PDEs. In this context, the 
development of the relevant apparatus for the asymptotic analysis of the oscillatory matrix RHPs 
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was originated in 1973-1977 in the works - ||3B|. It was essentially completed (for a detailed 
historical review see [^) in 1993 in the paper |^ where a nonlinear analog of the classical steepest 
descent method for oscillatory Riemann-Hilbert problems was suggested. 

The Deift-Zhou nonlinear steepest descent method consists of three basic steps (see see 
also |3^ and ||3j 



1. A deformation of the original jump contour S to the steepest descent (with respect to the 
dispersion exponent D{X; x, t)) contours and asymptotic evaluation of the solution away from 
the corresponding saddle points. 

2. The use of the relevant Lax pair and certain model Riemann-Hilbert problems to construct 
a parametrix for the solution near the saddle points. 

3. Assembling the above pieces into a uniform asymptotic solution which makes it possible to 



justify the whole construction by standard estimates |^0[ of the theory of singular integral 
operators on the complicated contours. 

Each of these steps has its natural analog in the classical steepest descent method for oscillatory 
contour integrals yet exploits much more sophisticated technics and analysis. In contrast with the 
classical steepest descent method, which is used for asymptotic evaluation of oscillatory integral 
representations, the nonlinear steepest descent method deals with a special type of oscillatory 
(singular) integral equations. In fact, the RHP ( |1.10| - |LT2|) is equivalent to the following singular 
integral equation: 



dfx 



2m Js — A+ 

The nonlinear steepest descent method provides a regular way of finding the proper transformation 
(which is highly nontrivial and virtually impossible to be seen directly!) of the original singular 
equation ( |1.15| ) to an equivalent one with uniformly small kernel. 

Some of the principal ideas involved in the steps 1, 2 of the nonlinear steepest descent method 
(e.g. explicit solutions for the model Riemann-Hilbert problems) go back to the earlier works 
and These earlier versions of the Riemann-Hilbert approach had already been successfully 
exploited in the asymptotic analysis of the various temperature correlation functions (see [^, ||29(| , 
Wli). The use of the nonlinear steepest descent method increases considerably the power of the 

Pi 



original scheme of |^, |r8|. In partricular, in |^0[ the use of the nonlinear steepest descent method 
allowed to compute the long-time asymptotics of the autocorrelation function of the transverse 
Ising chain at the critical magnetic field for the first time at finite temperature. The method has 
also produced solutions for some long-standing problems in the theory of random matrices and 
orthogonal polynomials (see e.g. [^, |Q, and @]). 



The XXO magnet and impenetrable Bose gas are free fermionic models. The non-free fermionic 
case is much more complicated because of several reasons. First, the Riemann-Hilbert problems. 
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which describe the corresponding Fredholm determinants, become operator- valued (cf. [^): 
the matrix elements Gj^kW turn into integral operators Gj^kW acting in an auxilary L2 space 
{Gj^kW Gj^fc(A|M, f )). In other words, an infinite-dimensional environment, absent in the 
free fermionic problems, arises. This transforms the associated Lax pairs and classical nonlinear 
PDEs into their non-Abelian analogies and obviously provides new significant difficulties for the 
asymptotic analysis. In fact, up to now the only attempt to solve an operator- valued RHP, which 



is related to correlation function of XX Z magnet, had been done in [^ |. 

The second difficulty is more subtle than the first one, and it is related to the presence of another 
infinite-dimensional context in the non-free fermionic problems; this time - of the quantum field 
nature. The fact of the matter is that out off free fermionic point the correlation functions can 
not be presented directly in the determinant form. Due to existence of non-trivial S-matrix, one 
needs to introduce auxiliary quantum operators — Korepin's dual fields — in order to find 
such a representation. As a result, the determinants obtained depend on these operators, and the 
correlation functions are equal to the vacuum expectation values of the Fredholm determinants 
in an auxiliary Fock space. In particular, for the correlation function ( p. . 1|) the following equation 
takes place (cf. [^), 

(*(0, Q)¥{x, t))T = const ■ e-^^'iOlBhtP, (j)A, (Pd],x, t) |0) , (1.16) 



where the factor const only depends on the temperature, coupling constant and chemical potential, 
and B is an operator acting in the auxilary Fock space with the vacuum vector |0). It involves 
the Fredholm determinant which functionally depends on the three basic quantum operators (dual 
fields) ■?/'(A), 0yi(A), and 0d(A). The exact definitions of the quantum operator B is given in 
appendix A. The dual fields ip^X), (PaW, (poi.^) are defined in section ^ by equations ( p.22| ) - 

dH). 

The operator- valued RHPs allow to evaluate the Fredholm determinants, but not their vacuum 
mean values. The two infinite-dimensional contexts, i.e. the operator nature of the RHP and the 
presence of the dual quantum fields, are completely unrelated. Therefore, in order to find the large 
time and long distance asymptotics of a non-free fermionic correlation function one needs to be 
sure that the asymptotics of the mean value is equal to the mean value of the asymptotics. 

It had been shown in Q that a naive asymptotic analysis of the Fredholm determinants con- 
taining dual fields does not provide a satisfactory result, i.e. the asymptotics is not uniform with 
respest to the avaraging over the dual fields. Thus, the problem arises — to obtain the asymptotic 
description of the determinant related to ( |1.1| ) which would be stable with respect to the procedure 
of averaging over the dual fields. This is the problem which we deal with in the present paper. 

Let us describe now the content of the paper. 

In section ^, following ||l| - 0, we give the basic formula and definitions concerning the non-free 
fermionic version of the determinant representation - Riemann-Hilbert approach to the asymptotic 
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analysis of the correlation function ( |1.1| ). 

In section ^, following we formulate the central object of the analysis, i.e. the operator- 
valued RHP associated with the correlation function ( p..l| ) (see ( p.lD and ( |3.11| ) below). The 
jump contour S of the RHP coincides with the real line, and the corresponding jump operator 
G{X) is realized as a 2 x 2 matrix whose entries Gjk{X) are the integral operators in an auxilary 
L2(— oo, oo) space. A principal feature of this RHP is that the operators Gjk{X) have the following 
special structure: 



^.^(A) = i6jk + \j){Gjk - Sjk){kl J, k = l,2, (1.17) 

where i is the identical operator in L2(— oo, oo), i.e. its kernel is a delta - function: i{u, v) = 6{u—v). 
The symbols \j) = \j,u) and {k\ = {k,v\ denote certain elements of L2(— oo,oo) (see definition 
( |3.9| ) below) satisfying the normalization condition, 

/oo rOG 
{l,u\l,u) du = {2\2) = {2,u\2,u)du=l. (1.18) 
-oo J —oo 

The numerical matrix G{X) in (|1 . 1 7|) is closely related to the jump matrix of the free fermionic 
impenetrable Bose gas model (cf. p9|). 

Equations similar to ( |1 . 1 7] ) take place for other non-free fermionic correlation functions as well. 



and they have a very important algebraic meaning. In fact, as it can be easiely verified, formulae 
(|1.17|) and ( |1.18|) define a representation. 



G^G = AiG), (1.19) 

of the group GL(2, C) in the group of Fredholm invertible operators in L2((— oo, oo), C^). It also 
can be shown (see |^^) that representation ( |1.19| ) preserves determinants, i.e. 

detG = detG. (1.20) 

The link ( |1.19| ) between the non-free fermionic jump operators and their free fermionic counterparts 
was first noticed by Korepin, and it has already proved very usefull in the asymptotic analysis of 
the non-free fermionic correlators (see |^ and |^). The mapping (|1.19|) plays a crucial role in 



the following sections where we develop the operator-valued version of the Deift-Zhou nonlinear 
steepest descent method. 

The main obstacle in taking the full advantage of the group-representation nature of equation 
(|1.17|) is that the vectors | j) and {k\ depend on the parameter A, i.e. | j) = \j, u, A) and {k\ = {k, v, A| 
(see ( p.9| ) below); moreover, they become singular for complex A. Therefore, one can not solve 
a non-free fermionic (operator ) Riemann-Hilbert problem by just applying the representation 
(|1.19|) to the solution of the corresponding free fermionic (matrix) Riemann Hilbert problem. 
Nevertheless, this representation helps to perform the first step of the nonlinear steepest descent 
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method (the contour deformation) and evaluate the leading term of the asymptotic solution of the 
RHP. This is done in sections ^, ^ where we closely follow the methodology of . The leading 



asymptotic term we found coincide with the one obtained earlier in via the direct asymptotic 
analysis of the Fredholm determinant. 

The objective of sections ^j-^ is an order t^^/^ correction. This is the second step of the 



nonlinear steepest descent method. In the free-fermionic case (see e.g. ||3^ and ||3^) this step 
constitutes the reduction of the deformed RHP to a model problem associated with the saddle 
point. The model problem is then solved explicitly in terms of the parabolic cylinder functions (see 
agam m and [|^; see also H). In our case, the model problem is an operator-valued version of 



the classical model problem. In fact, the jump operator of the model problem is the image of a 
certain jump matrix, closely related to the corresponding free fermionic model jump matrix, under 
the modified mapping Aq. The latter is the representation A with the vectors \ j) and {k\ replaced 
by their values evaluated at the saddle point Aq, 

\ j,u) = \j,u,Xo), {k,v\ = {k,v,Xo\. (1.21) 

The representation Aq does not depend on A and hence can be used to transform a free fermionic 
matrix model solution into the operator model solution. The transfortmation though is not quite 
straightforward, but it allows eventually (see sections |^, H) to solve the model RHP explicitly 
(more exactly, up to the inversion of an integral operator which is independent on A, x, and t). 
In achieving this result, the following 'operator - indexed' generalization Dc,{C) of the classical 
parabolic cylinder functions D^{$,), 

DoiO = DoiO - I hO l) + I J){J I MO, ^ = ^1 J){J I , (1-22) 

plays an important role. 

The order correction is a main threshold in the asymptotic analysis of the temperature 
correlation functions (cf. e.g. |^). After it is passed, one can, in principle, obtain a total 



asymptotic expansion using the elementary algebra only, just substituting a proper asymptotic 
series into the nonlinear PDEs associated with the correlation function under consideration. In 
our case, the relevant nonlinear PDE is the mentioned above non-Abelian version of the classical 
NLS equations (see (|2.21| ) below). In section ^ we use this non-Abelian NLS for exact evaluation 



of a few terms next to the order t^^/^ term and for evalution of the pre-exponential factor in 
the asymptotics of the Fredholm determinant. The non-Abelian NLS also allow us to extract a 
valuable information about the full asymptotic expansion. We analyse this information in section 
and conclude that to make the asymptotic expansion obtained stable with respect to the dual 
field avaraging a certain modification of our construction is needed. The modification concerns 
the equation determing the saddle point Aq, and it is made in the last section |T^. 
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As the main result of the paper, we suggest the following asymptotic formulae for the quantum 
operator B in the r.h.s. of equation (|1.16|) , 



i3 = C_(0^,0A|Ao,T, h,c){2ty 



X 



exp /°° - 2At + sign(A - A) 

: ln{l - i9(A) (^1 + e^W -s-(^-a)^ } dxj (l + O (^-^ 



t —y oo, 



X 

2t 



= Xo = 0{l] 



for the case of negative chemical potential h, and 



;i.23) 



B = 0^1 Ao, T, h, c)(2t)-— e 



^^pV'(Ai)+itAi2-ixAi 



exp 1^ ^ (^x - 2 At + ttp'iX)^ sign (A - 3?A) 



XA) sign{KA-A) 



t^oo, - = Ao = 0(l), 



:i.24) 



for the case of positive chemical potential h. 

A few remarks explaining the notations which are used in equations (|1.23|) , ( |1.24|) and the 
meaning of the equations themselves are needed: 

1. The function 'i9(A) is the Fermi weight ( |l.(j| ). 

2. The dual field 0(A) equals the difference 0yi(A) — (pn^X). 

3. The exponent i^(A) is defined by the formula, 

1 



27Ti 



In 



i9(A)(l + e-^(^))) (l - ^9(A)(1 + e'^(^))) 



4. Ai and A are the roots of the equations 

l-^9(Ai)(l + e-'^(^i)) 



0, /i>0, 



:i.25) 
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and 

A = Ao + ^^'(A), (1.26) 
respectively. Since A is not necessary real one should understand the integrals ( |1.23| ), ( |1.24| ) 



as 



/OO P .A. /* OO 

F(sign(A-A),A)rfA= / F{l,X)dX+ F{-l,X)dX. (1.27) 
-OO ^ ^ J—oo J A 

5. The undetermined factors C±{(j)A, 4>d\^o,T, h, c) are functional of the dual fields (pA, 4>d, 
and they are functions of the temperature T, the chemical potential h, and the coupling 
constant c. They do not depend on the dual field ip, and they only depend on distance x 
and time t through the ratio Xq = x/2t. 

6. All the commutation relations involving the dual fields 0a (A), ^^(A), and ipW ^i^re trivial (see 
( 2.24 )). Therefore, when deriving ( |1.23 ) and ( 1.24| ) we treat the dual fields as the complex- 



valued functions analytic in the strip |ImA| < c/2; the functions exp0A,D and exp-^ are 
supposed to be bounded in the strip. It is also assumed that equation ( p..25| ) has no real 
solutions a h < 0, and it has exactly two real roots, i.e. Ai, A2, if h > 0. Moreover, it is 
supposed that 

Ai < A2 < A, (1.28) 

and that the same inequalities, in the case h > 0, take place for the roots of the function 
1 — 'i?(A)(l + e'^*^'^^). These assumptions are justified in sections 0, ^ and ^ Accordingly, the 



contour of integration T in (|1.24|) is choosen as it is indicated in Figure 2 below, so that the 



integral in ( p.. 24 ) makes sense. 



7. The point A, which plays the role of the 'shifted' saddle point, is defined by equation (|1.26|) 
as its root which approches Aq as t ^ 00. In other words, A — Aq = 0{t~^). Nevertheless, 
we do not replace A by Aq in the asymptotic formulae ( |1.2^j| - 1.24| ). As it is explained in 



sections ^ and |TT|, to make the asymptotics of B compatible with the avaraging over the dual 
fields we need to trace the presence of the field in all the corrections to the leading terms 
indicated in ( p..23| - [1.241 ). In fact, we need all the corrections to be of the order 0{ip'^/t'^) 



with < m < n. The replacement A Aq, if made, would introduce the corrections of the 
order 0{ip™' /t^) with m > n which produce the positive powers of t after the avaraging in 
(|1.16| ) (see sections |10| and |Tl| for more details). In other words, when the operator nature 



of equation ( |1.26|) is taking into account one has to keep in mind that, although the factor 



i/2t is small, the quantity ip'i^X) is an unbounded operator of the derivative type. 

8. Because of the quantum operator nature of the dual fields (Pa{X), 0d(A), and ip{X), the r.h.s. 
of equations ( |1.23| ) and ( |1.24| ) have, strictly speaking, a symbolic meaning. One needs the 



method of avaraging of the functionals of the dual fields developed in [B in order to evaluate 
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the mean values of the objects indicated in ( |1.2ij| - |1.24D . The results of these calculations, 
which would finalize the asymptotic analysis of the leading terms of the correlation function 
( |l.lj ), will be given in the forthcoming paper. 



We conclude the introduction by emphasizing that in this paper we do not address the questions 
related to the third step of the nonhnear steepest descend method, i.e. the questions related to the 
rigorous justification of the asymptotic equations ( [1.23| - |1.24| ). Some of these questions, including 
the rigorous setting and the general theory of the operator-valued Riemann- Hilbert problems 
arising in the non-free fermionic models, are still under the investigation. The purpose of this 
paper is to show that, using the operator- valued Riemann-Hilbert problems, it is posible to carry 
out an effective asymptotic analysis of the general time-dependent temperature non-free fermionic 
correlation functions based on their determinant representation. 



A = A{u,v) 



(2.1) 



2 The basic formulae and definitions 

This sections contains the basic definitions and notations, which are used in the paper. The reader 
may find the details in the papers [jT|, 0. Here we restrict our selves with only necessary formulae. 
Consider the operators of the following form 

^ An{u,v) Au{u,v) 

\ A2l{u,v) A22{U,V) 

The entries of this 2x2 matrix are integral kernels Ajk{u,v) acting on functions from the 
L2(— oo,oo) space. We shall mark these operators (as well as their entries) with the symbol 
'hat'. Also, we will frequently use the equation, 

M = M{u,v), 

to indicate that M{u, v) is the kernel of an operator M. 
The product of two operators of this type is defined as 

iAB)jkiu, v) = Y, Ajiiu, w)Bikiw, v) dw. (2.2) 

We shall also consider the products of separate matrix elements: 



oo 



{AjkBim)iu,v) = Ajk{u,w)Bi^{w,v)dw. (2.3) 
The definition of trace of 'hat '-operators is standard 

Tri = trill +tri22 = / (Au{u,u) + A22{u,u)) du, (2.4) 

J — oo ^ ^ 
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where 

tr Ajk = I Ajk{u,u)du. (2.5) 



The main object of study in this paper is an operator-valued Riemann-Hilbert problem. Here 
we give the general formulation of this problem. 

The RHP consists of the following: find an operator xW depending on complex parameter A 
and satisfying the following conditions (cf. - ( 1.12|) ): 



x(A) ^ /, A ^ oo, (normalization condition) (2.6) 
2^°\ X{X) is analytical function of A if A ^ i?, (2.7) 
3(0). x-(A) = x+(A)G'(A), XeR, (jump condition). (2.8) 



Here / is the identity operator 




6{u-v). (2.9) 



Operators x±{X) are boundary values of the function x(A) from the left (from up half-plane) and 
from the right (from low half plane) of the real axis respectively (compare again with (|1.1CI| ) - 
(|1.12|) ). G'(A) is a given operator. A particular choice of this operator which corresponds to the 
correlation function ( p. . 1| ) will be described in detail in the next section. 

Following the tradition we call the operator G{X) 'jump matrix', although this matrix evidently 
is infinite dimensional. We hope that this terminology will not cause a confusion. 

Although the RHP (0) looks like 2x2 matrix RHP, here we deal with infinite-dimensional 
integral operators. In more detail, the jump condition 3^^^ can be written as 

^ poo ^ 

Xjk-{X\u,v) =J2 Xji+{X\u,w)Gik{X\w,v)dw, XeR. (2.10) 

The analyticity of the operator x(A) is understood as the A-analyticity of its kernel x{X\u,v). All 
the limits and expansions involving x{X) = x{X\u,v) are supposed to be uniform with respect 
to u and v. We will not go further in the specification of the setting of the RHP (|2.6| ) - ( p^ 



As it has already been indicated in the introduction, in this paper we do not discuss the general 
theory of the operator-valued RHPs. In what follows, we assume that the unique solutions of 
the operator-valued RHPs we are dealing with always exist and satisfy all the natural properties. 
For instance, we shall assume that if G{X) is a Fredholm operator, i.e. the determinant det(j'(A) 
exists, then x{X) is a Fredholm operator, and its determinant detx(A) is a unique solution of the 
following scalar RHP, 

1. detx(A) ^1, X^oo, (2.11) 

2. detx(A) is analytical function of A if A ^ i?, (2-12) 
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3. detx-(A) = detx+(A)detG(A), X e R. 



(2.13) 



In particular, this means that the function det x{X) is given by the exphcit formula, 

detx{X) = expi^-^J"^^-^\ndetG{fi)'^ , (2.14) 

subject the zero-index solvability condition, 

argdetG(A) ^ 0, A ^ ±oo . (2.15) 

The asymptotic expansion of the solution of the RHP at A ^ oo will play the central role 



x(A) = /+^ + ^ + ... . (2.16) 



The first two coefficients of this expansion b and c allow one to reconstruct the Fredholm determi- 
nant, i.e the r-function corresponding to the solution xW: which in turn directly related to the 
correlation function of the local fields ( |1 . 1|) . This correlation function can be presented as 

(^(0, 0)¥{x, t))T = const ■ e~"'\0\B{[tP, 0^, <Pn],x, t) |0). (2.17) 

We have extracted in this formula the factor exp{— i/it}, which trivially depends on t, and a 
constant, which does not depend on x and t. Strictly speaking this 'constant' depends on the 
temperature, coupling constant and chemical potential. The most important object in ( 2.17|) 



is By[ip,(j)A,(f)D],x,tj, which, in particular, includes a Fredholm determinant, and which is an 
operator in auxiliary Fock space. This operator is a function of x and t and it functionally 
depends on three quantum operators (dual fields) i^^X), 4>a{X) and 0d(A). The last ones act in 
the auxiliary Fock space having the vacuum vector |0) and dual vector (0| (one should not confuse 
the auxiliary vacuum vector |0) with the vector |0) in ( |1.3| )). Thus, the correlation function is 
proportional to the vacuum expectation value of the operator B. The detailed definition and the 
properties of the dual fields will be given later in this section. First we want to focus our attention 
on the relationship between the operator B and the operator-valued RHP. 
The following representation for the operator B had been found in [|ll, P 



B 



t) = det {l + v) ■ p bu{u, v) dudv. (2.18) 



We call ( p.l8| ) 'determinant representation', since the r.h.s. of this formula is proportional to 
the Fredholm determinant of the linear integral operator I + V. The detailed description of this 
integral operator is given in Appendix A. 
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The results of the paper |^ allow one to express the Fredholm determinant and factor 612 (m, v) 
in terms of the solution of an operator- valued RHP (0) with a certain jump matrix G{X). In fact, 
the operator 612 in ( p.l8| ) is the corresponding entry of the operator coefficient b of the asymptotic 



expansion ( p.l6|) . The logarithmic derivatives of det(^/-|- with respect to distance x and time 
t also can be written in terms of the coefficients b and c from ( |2.16| ): 

9^.1ogdet(/ + =itrbi 



ni, 

dt log det i^I + Vj = iitr(c22 - cii). 
The second logarithmic derivatives depend on the 'matrix elements' 612 and 621 only: 
dAhgdet(i + V) = -tiibuki), 



(2.19) 



. . . . (2-20) 

dtd^ logdet^J + Vj = Ur(9^6i2 • 621 - <9x&2i ■ bu). 

Finally these two operators satisfy a non-Abelian generalization of the classical Nonlinear Schro- 
dinger equation 

-idtbu = -dlbi2 + 2612621&12, 

(2.21) 

idtb2i = -dlb2i + 2621612621- 



Recall that one should understand the products of operators 612 and 621 in the r.h.s. of ( f^.20D and 



(12:211) in the sense (|J) 

Thus, we have described E in terms of the solution of the operator-valued RHP (|27 
(the explicit expression for the jump matrix G{\) will be given in the next section). Equations 
( |2.19|) define ;B up to a constant factor. Equations ( |2.20|) , (p.21|) will help us to analyse the higher 
corrections to the leading asymptotics of 612 and 621. 

Let us turn now to the auxiliary quantum operators — dual fields. These operators ip{\), 0_d(A) 
and 0yi(A) (originally the last two fields were denoted as 0z)i(A) and (t>A2W) were introduced in 

in order to remove two-body scattering and to reduce the model to free fermionic one. As it 
was mentioned already, they act in an auxiliary Fock space. Each of these fields can be presented 
as a sum of creation and annihilation parts 

0a(A) = gA(A)+PD(A), 

0d(A) = gz5(A)+PA(A), (2.22) 
^(A) = g^(A)+p^(A). 

Here p(A) denotes the annihilation parts of the dual fields: p(A)|0) = 0; g(A) denotes the creation 
parts of the dual fields: (0|g(A) = 0. 
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The only nonzero commutation relations are 

[pA{X),qM] = = ln/i(/i,A), 

[Pd(A), = [p4X), qDifi)] = In h{X, (2.23) 

[Pi>{X), = ln[/i(A, fi)h{n, A)], where /i(A, /i) = - — ^ilLl£. 

Recall that c is the coupling constant in (p..5|). It follows immediately from (|2.23|) that the dual 
fields belong to an Abelian sub-algebra 

[V^(A),V^(/i)] = [V^(A),0„(/i)] = [MX),Mf^)] = 0, (2.24) 



where a,b = A,D. Due to the property (|2.24|) , the Fredholm determinant det(^J + V^^ is well- 
defined. 



The Fredholm determinant det(^/ + V^^ and the factor bu in the representation ( p.l8|) func 



tionally depend on the dual fields. Hence B is an operator in auxiliary Fock space (for explicit 
formulae see ( |A.1| )-( |X^ )). From the equations ( ^.191 ), ( ^.2(J| ) one can conclude, that the solution of 



the corresponding operator-valued RHP, as well as the jump matrix, also should depend on these 
auxiliary operators. We shall see in the next section that the jump matrix G{X) does depend on 
dual fields, therefore it is an operator in the auxilary Fock space as well. 

However, we would like to emphasize that the operator nature of the RHP (0) is not related 
to the fact that x{X) and G(A) depend on dual fields. Even if we replace all dual fields by some 
complex functions, the RHP (0) still remains the operator-valued one, since the entries of x{X) 
and G{X) are integral kernels. 



The calculation of the asymptotics of the correlation function ( |1 . 1| ) consists of two stages. At 
the first stage one has to solve the operator-valued RHP (0) in order to find explicit expression 
for the operator B. The second stage consists of the averaging with respect to auxiliary vacuum 
vectors. The role of dual fields is significantly different at these two stages. While the averaging is 
based essentially on the operator properties of dual fields, at the first stage the analytical properties 
of the last ones are much more important. While solving the RHP (0) one can consider the dual 
fields as some complex functions, which are holomorphic in a neighborhood of the real axis. Indeed, 
all these auxiliary operators commute with each other, and their matrix elements are holomorphic 
functions if their arguments belong to the strip | ImA| < c/2. 

However such a treatment of dual fields leads us to a problem, which we had touched briefly 
in the Introduction. The matter is that the operator-valued RHP considered in the present paper 
can be solved only asymptotically for the large time t and long distance x. Thus, we obtain only 
asymptotic solution for RHP, and, hence, we find only the asymptotics of the operator B. Strictly 
speaking, one has to be sure that vacuum expectation value of the asymptotics is equal to the 
asymptotics of the vacuum expectation value. Otherwise corrections to the asjTiiptotic expression 
for B may give non- vanishing contribution into (0|i3|0). 
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Below (section |TT|) we shall see that there exists a set of different asymptotic expansions of 
the operator B, which are equivalent, if we treat the dual fields as complex functions. All of 
them provide us with the same result for the vacuum expectation value if and only if we take into 
consideration the complete asymptotic series. However, if we restrict our selves with the leading 
term of asymptotics and finite set of corrections, then these asymptotic series give different results 
for (0|i3|0). In section [Tl| we will suggest an algorithm of choosing the asymptotic expansion which 
is compatible with the avaraging over the dual fields. 

In summary we can formulate our point of view on the dual field problem as follows. We 
consider the dual fields as the complex functions, which are holomorphic in a finite neighborhood 
of the real axis. Moreover, since the matrix elements of exp0a(A) and expiplX) are rational 
functions of A, we shall consider exp(/)a(A) and exptp{X) as the rational functions of A, bounded 
in the strip | ImA| < c/2. This permits us to solve the operator-valued RHP for the large time 
and long distance separation by deforming the original contour to the contours in the complex 
plane. At the same time we do not forget about the quantum operator nature of these auxiliary 
operators. Therefore we do not restrict our selves with the leading term of the asymptotics for the 
operator B. We draw our special attention to the corrections and their behavior after averaging 
with respect to the auxiliary vacuum. We shall consider all these questions in more detail in 
sections 10, 11, where we will demonstrate how one can modify the procedure of solving the RHP 
in order to obtain correct result for the asymptotics of the vacuum mean value. 



3 The operator-valued RHP 

The operator-valued RHP describing the temperature correlation function of local fields (|1.1|) was 
formulated in [Q. We have presented already the general form of this RHP in the previous section. 
Here we give the explicit expression for the jump matrix and explain all necessary definitions and 
notations. 

According to 0, the operator- valued RHP in question consists of finding the operator x{X) 
possessing the properties, 

r. x(A) ^ /, X^oo, 

1°" . x(A) is analytical function of A if A ^ i?, (3.1) 
3^X-(A) = x+(A)G(A), AGi?, 
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where the operator-valued entries of the jump matrix G{X) are defined by the equations, 

GuiMu,v) = 6{u-v)-Z{v,X)6{u- X)^{X)e'^^^^^] 

GuiMu, v) = 2m{^{X) - l)6{u - X)6{v - A)e^(^)+"(^); 

7 (3.2) 

G2i{X\u,v) = --^Z(M,A)Z(t;,A)?9(A)e'^^(^)+'^«(^)-^(^)-"(^^ 

ZTT 

G22iX\u,v) = 6{u-v)-Z{u,X)6{v-X)^{X)e''''''^^\ 



Here (Pa{X), 4>d{X) and i'iX) are the dual fields ( p.22| ). Recall that starting from this section 



we consider exp0^(A), exp0£)(A), and expip{X) as the classical functions which are rational and 
bounded in a finite strip ( |ImA| < c/2) near the real axis. Moreover, we will assume (see the 
arguments given in the end of this section ) the following symmetric properties of the functions 
7/>(A) and 0(A) =0^(A)-0z,(A) : 

^(A)=^(A), (3.3) 

and 

4>{X) = -0(A), (3.4) 

The function Z{X,fi) is equal to 

p-<t>DW p-<I>aW 

where h{X, fi) = {X — fi + ic)/ic was introduced in ( |2.23| ). The Fermi weight '(9(A) was defined in 
(|1.8|) , and it is rapidly decreasing function at A — > oo. The only object depending on the time t 
and the distance x is the function r(A) = r(A; x,t) defined by the equation : 

r(A) = itX"^ - ixX, (3.6) 

and representing the dispersion law of the underlying classical model. 

In the present paper we study the large time and the long distance asymptotic behavior of the 
solution of the RHP (a). More precisely we consider the case when t oo, x ^ oo, while their 
ratio remains fixed: 

x/2t = Xo = 0(1). 

Therefore it is convenient to use new independent variables t and Aq instead of t and x. The 
function r(A) then turns into 

r(A) = it{X - Ao)^ - itXl (3.7) 

Since the jump matrix G{X) ( p.2|) depends on delta-functions whose arguments contain pa- 
rameter A, one has to understand the jump condition 3" in a weak sense. In order to avoid the 
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complications caused by the presence of the distributions depending on the complex parameter, 
we make the following regularization of the jump matrix ( p.2|) . 
Let 

5,{u - A) 



1 (u-A)^ 



2^/7^e 



(3. 



and 



M(A) 



5^{u — \)Z{u, A) 
'M(A)Z(A,A) ^ 



A)(5e(M — A) (i-u, 



Z(A,A) 



\2) = \2,u) 



Z{\\) 



ZKA), (2| = (2,t;| 



\ M(A) 

A)Z(i;,A) 
y'A4(A)Z(A,A) ■ 



(3.9) 



Obviously 



\ M(A) 

/OO /"OO 
(l,M|l,M)rfM= (2|2) = / (2,u|2,u)du 
-oo J —oo 



(3.10) 



one can see that AC (A) has no zeros on the real 



Also, in virtue of the symmetry equation 
line. 

We define the entries of the regularized jump matrix by the equations, 

Gn(A) = ?-^(A)Z(A,A)e'^-(^)|l)(l|; 

G'i2(A) = 27r2(^9(A) - l)Z(A,A)e^(^)+"(^)|l)(2|; 

G'2i(A) = -— 7?(A)Z(A,A)e'^^(^)+^^(^)-^(^)~"(^)|2)(l| 
27r 

G22(A) = ?-^(A)Z(A,A)e*^W|2)(2|, 



(3.11) 



5{u 



It is easy to see that in the limit e ^ the jump matrix (^.11|) turns into 



where i 

dU). 

The RHP (a) with the regularized jump matrix corresponds to a new regularized operator 
B. In particular, the Fredholm determinant + should be replaced by + (see 

appendix B). New regularized kernel V^, as well as original kernel V , is not singular, and the limit 
— > e — > is well defined. It is remarkable that all relationships (|2.19|) -( ^.21|) remain valid 
for the regularized objects for finite values of e (not only in the limit e — > !). In the rest of the 
paper we will only work with the regularized operator V^, and we will denote it as V omitting the 
subscript e. We will also assume that the large time limit we study commutes with the limit e — > 0. 
As an indirect justification of this conjecture, we take another special property of our regularization 
which will be revealed later: the first leading terms which we evalute for the regularized Fredholm 
determinant do not depend on e. 



18 



As in under the assumptions made in the previous section concerning the general theory of 
the operator- valued RHPs we are dealing with, solvability of the RHP (|3.1| , p.ll|) implies uniqness 
of its solution. Moreover, using the representation property ( p..2q ) and equation ( |2.14| ) we conclude 
that the solution xW of the RHP (a) with the jump matrix (|3.11|) satisfies the equation. 



det x{\) = 1 for all A . 



(3.12) 



The very possibility of solving the RHP ( |3.1| , |3.11|) asymptoticly as t ^ oo is based on the 



presence of the rapidly oscillating exponents e^^^'^-' in the jump matrix. In the next section we 
start the asymptotic analysis of the solution of the RHP { ^A\ , p.ll|) using an operator-valued 



generalization of the nonlinear steepest descend method suggested in for the oscillatory matrix- 
valued Riemann-Hilbert problems. 

We conclude this section by one important note. 

The model QNLS has two different phases, corresponding to the positive and negative chemical 
potential h in the Hamiltonian ( p..5|) . In particular, the ground state of the model for < 
coincides with the bare Fock vacuum, while for > the ground state is the Dirac sea At 
finite temperature the difference is not so essential, but nevertheless one could expect that the 
asymptotics of the correlation function is different for these two cases. Since this asymptotics is 
defined by the solution of the RHP (a), the properties of the jump matrix ( |3.11| ) should be different 



for h > and h < 0. This fact does take place in the free fermionic limit and it is related 
to the zeros of the diagonal elements of the jump matrix. Out off free fermionic point one should 
consider zeros of the determinants of Gu and G22- Since both of these operators are the sum of 
identity operator and one- dimensional projector, it is easy to see that 



det Gii(A) = 1 - ^iX)Z{X, A)e^^(^) 



e T -\- I 

(3.13) 

det G22(A) = 1 - mZ{X, A)e^-^(^) = —^^ , 

e T + 1 

In the free fermionic point (coupling constant c goes to infinity) one can put dual fields in these 
formulae equal to zero, since due to commutation relations ( p.23| ), the vacuum expectation value 
is trivial in this limit. As we had mentioned in the Introduction, the solution of the Yang- Yang 
equation ( |1.6| ) £:(A) is positive, if /i < 0, and £:(A) has two real roots, if /i > 0. Therefore the 
determinants of Gu and G22 possess just the same properties. 

Out off free fermionic point, in order to answer the question, whether the determinants of Gu 
and G22 have zeros at the real axis, one need to know the properties of the function exp{£:(A)/T} — 
exp{±0(A)}. These properties must reflect by a natural way the properties of the operator 
exp{e(A)/T} — exp{±0(A)}. It was shown in [1] that dual fields can be expressed in terms of 
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canonical Bose fields. Using these representations one can prove, in particular, that V^^(A) = V'(A) 
and 0^(A) = —0(A). Therefore we demand the corresponding classical functions to possess similar 
properties ( p.3| ), (p.4|). The operator exp{±0(A)} obviously is a unitary operator with the spec- 
trum, belonging to the unite circle. There are serious reasons to assume that all the points of the 
unite circle belong to the spectrum of this operator (although this fact have not been checked rig- 
orously). Therefore, if h < 0, and hence e(A) > 0, we conclude that the spectrum of the operator 
exp{£:(A)/T} — exp{±0^(A)} does not contain zero: 

^ CT (e^(^)/^ - e±<^(^)) , h<0. (3.14) 
On the contrary for h > the function e{X) has two real roots, hence 

^ a (e^^yT _ e±m) ^ /,>0. (3.15) 



Thus, taking into account ( |3.14| ), ( p.l5[ ), we shall demand the classical function 0(A) = 0a (A) — 



0£)(A) to be such, that the determinants of the operators Gu and G22 have no real roots for h < 0, 
and they both have two real roots for /i > 0. 



4 Negative chemical potential 

In this section we begin the asymptotic analysis of the RHP (a) with the jump matrix ( ^.llj ). We 
refer the reader to the works ||3^ and which present the nonlinear steepest descent method for 
the classical NLS, to see that each of our basic operator constructions has its matrix counterpart. 
In particular, in this and the next sections we perform the first step of the method, i.e. the 
transformation of the jump matrix to the proper up-and-low-triangular forms followed by the 
deformation of the real line to the steepest descent contours with respect of the exponent t(A). 



All the considerations of this section are very close to those of |^ . 

We start with the case of negative chemical potential. Let as before Aq denote the saddle point 
of the exponent t(A), i.e. (see (|3.6| ) - (|3.7|) ), 



A„ = |. (4.1) 

and consider the following substitution 

X(A) = <l(A)7t(A), (4.2) 
where $(A) is new unknown operator- valued matrix, and 7?.(A) is diagonal matrix 

m 

n=\ I . (4.3) 

(^-(A))- 
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Here ^(A) is the solution of the RHP 

1^. ^(A) — > ? = 5{u — v), A — > oo, 

2^. ^(A) is analytical function of A if A ^ i?, (4.4) 

3\ g^X) = g+{X) (9{\o - X)Gu + ^(A - Ao)(G'^2)"') , A G i?, 



and its kernel g{X\u,v) is integrable at the point A = Aq (cf. the function (5(A) in [^], [p7| ). 
The RHP (6) looks like a scalar RHP, however one should remember that ^(A) = g{X\u,v) and 
Gjk{X) = Gjk{X\u,v) are integral operators. The symbols ^"^(A) as well as G22W mean the 
transposition of these operators: g^{X\u,v) = g{X\v,u) and G22iM'^y'^) = G22{X\v,u) respectively. 
9{X) is the step function 

dW = I ^' ! ^ °' (4.5) 
^ ^ ^ 0, A<0. ^ ^ 



Since the jump operator in ( [4.4|) is discontinuous at A = Aq we have demanded an additional 
condition of integrability for the solution of the RHP (6) (cf. the standard conditions fi^] of the 
theory of factorization of matrix functions). Recall also, that we consider the case, when det Gn 
and det 622 are not equal to zero at the real axis, hence the jump operator ^^(Ao — A)Gii(A) + 6{X — 
Xo)(G22{X)^ is well defined. 

Thus we obtain a new RHP for $(A) 

r. <|)(A)^/, A ^00, 

2"^. <I>(A) is analytical function of A if A ^ i?, (4.6) 
3^ <|)_(A) = <I+(A)G$(A), XeR, 

where 

G^{X) = 7^+(A)G'(A)7^:^(A). (4.7) 

The coefficients b and c of the asymptotic expansion of the original operator x(A) (the RHP (a)) 
can be easily expressed in terms of asymptotic expansions of $(A) and ^(A): 



g = t+j + j^ + .... (4.9) 



In particular, it is easy to see that 
hi = ($0)11 + ^0, 



2 (4-10) 

C22 - Cii = ($1)22 - ("^'Oii - ("^'0)22^^ - ('^'0)11^0 - Qi- + {qo) ■ 
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Consider the function A(A) = det ^(A). This is not an operator but a complex function. The 
asymptotic expansion of this function at A — oo plays an important role. It is convinient to write 
down this expansion in the exponential form, 



A = exp^^ + ^ + ,,,^ (4.11) 



In order to find the first logarithmic derivatives of the Fredholm determinant we need to know 
only the traces of the operators bu and C22 — Cn (see ( p.l9|) ). Therefore taking the trace of ( |4.10|) 
we arrive at 



where 



and 



where 



(4.12) 



tibn =tr($o)ii + Ao, 

tr(c22 - cii) = tr [($1)22 - ($1)11 - ($0)22^0 - ("^0)11^0] - 2Ai. 
Consider now the jump matrix Using the equations, 

(|l)(l|f = |2)(2|, {\k){j\r=\k}{jl k^j, (4.13) 

and 

m{m\j){i\) = \k){i\, k,j,i = 1,2, (4.14) 

one can easily check that 

G22 — (521 (Gil) G12 = (Gil) , 

Gil — Gi2{G22) ^G2l = {G22) ^■ 

Due to these properties the jump matrix (5$ has the form 



(4.15) 



G^=\ I , for A < Ao (4.16) 

Qe-' i + Qp 



Qe~^ = if+r'G2i{Gii)-'g-\ (4.17) 
Pe^ = ^-(Gii)-^Gi2^^, (4.18) 



i + PQ Pe^ 

G$ = I . I , for A > Ao (4.19) 



Qe-^ = if_r\G22)-'G2igZ\ (4.20) 
Pe^ = Q+Gi2{G22)-'gl. (4.21) 
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= $ 




u = <^> 



Figure 1: Contour for new RH problem 



Therefore one can factorize the jump matrix 

G$ = M+M_ 



where 



A^4 




for A < Ao, 
for A > Ao, 



M_ 



Qe 



(4.22) 



(4.23) 



(4.24) 



Notice that Q and P can be analytically continued into a neighborhood of the real axis in the 
upper half-plane. Similarly, Q and P can be continued into a strip in the lower half-plane. 

Consider new contour, which is shown on the Fig.[^ The contour consists of four branches Ci, 
6*2, C3 and C4, having the origin in the saddle point Aq (arrows show the positive direction). The 



branches Cj and the real axis are the boundaries of the four half- infinite 'wedges' Qj, {j 



.4). 



Matrices N± can be analytically continued into the wedges Qi and ^4 respectively, matrices M± 
— into the wedges fl2 and fl^. 
Define the matrix If as follows: 

f/ = <l>, A ^ u u u 1^4, 



U = ^N+, 
U = <I>M+, 

U = ^MZ\ 
U = ^N-\ 



A G fii, 
A e 1^2, 

A G fis, 

A G 



(4.25) 
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(see Fig.|Tp. Then the matrix U has no jump at the real axis, however it has jumps at the new 
contour C = Ci U C2 U C3 U C4. The RHP for U can be formulated now as 



l*^. f/(A) ^ J A ^ 00, 

2*^. U (A) is analytical function of A if A ^ C, 





= u, 




AgCi, 




= u, 




AeC2, 


u_ 


= u, 


.M_, 


AGC3, 




= u, 




A e C4. 



(4.26) 
(4.27) 
(4.28) 
(4.29) 



where are boundary values of If from the right and left of the contour C respectively. 
It is easy to see that all jump matrices in ( [4.26| )-( |4.29| ) have the form 

/ + o(e-*'*), for lA-Aol >ri/'+^ 



(4.30) 



where 6 is an arbitrary small positive number. Thus, we can put N± = M± ^ J up to exponentially 
decreasing corrections for all A G C except small vicinity of the saddle point. Formally appealing 
to the integral equation ( |1.15| ), written for the contour C and the operator- valued function f/+ (cf. 



the nonformal arguments of for the usual matrix case), we conclude that 

f/ ^ /, ^ !> ^ /, ^ X ~ 7^. (4.31) 
and hence the coefficients of the asymptotic expansion $0 and $1 are asymptoticly equal to zero: 

$0, $1 ~ 0. (4.32) 



Thus, similar to |^, the traces of the operators hu and C22 — Cn asymptoticly depend only on 
the coefficients of the A - expansion ([4.11|) of the function A(A) (see (|4.12| )): 



(4.33) 



trfen = Ao + o(l), 
tr(c22 -cii) = -2Ai + 0(1). 
On the other hand, one can easily find A (A) since the RHP (6) for the operator ^(A) implies the 



following scalar RHP for its determinant A(A) (cf. ( p.ll|) - ( |2.13|) ): 

r. A(A)^1, A^cx), 

2^. A(A) is analytical function of A if A ^ i?, 

3^ A_(A) = A+(A) det (^(Ao - A)Gii + e{\ - Ao)(G: 
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(4.34) 



T \-l 
22) 



Using the fact that both of the operators Gn and G22 are a sum of the identity operator i and an 
one-dimensional projector we find 



and (cf. ( CT ) 

det (^^(Ao - A)Gn(A) + ^(A - Ao)(Gi;(A))"') = (l - ^(A)(l + e<^W^'^°(^-^o)))^'s°(^-^) ^^ ^^^ 



where sign(A) = 6{\) — 6'(— A) is the sign function, and 0(A) = 0a(A) — 0d(A). Thus from ( |2.14|) 
it follows that the solution of the RHP (e) is given by the explicit formula 

A = exp|--^ ^^sign(Ao-/i)ln(l-i9(^)(l+e'^(^)^'S"(^-^°)))|, (4.37) 



2TTi J-oo /i — A 

and, hence, 

Ao = — / rf/isign(Ao - n) In (l - ?9(/i)(l + e<^(^)«is>^(^'-^'') ) ) , (4.38) 
27ri J-oo ^ ^ 

1 /"OO , , , ^ 

Ai = — / /irf/i sign(Ao - /i) In (1 - ^9(/i) (1 + e^^'^) sign(M-Ao)A \ ^ ^4 39) 

27r2 J-OO ^ ^ ' ' 

The logarithmic derivatives of the Fredholm determinant with respect to the variables t and Aq 
are equal to 

^aAolndet(/ + \/) =ztr6n =zAo + o(l), 

/ A \ ^'-'^^ 
[dt - y^Ao j lndet(/ + l^) = 2tr(c22 - hi) = -2zAi + o(l). 

Integrating these equations with respect to t and Aq we arrive at 

lndet(/+r) = — / dfx\x-2fxt\\n(l-^{fi)(l + e'^^''^''^''^'''^°n) +0{\nt), (4.41) 

2n J-oo ^ ^ 

where x = 2tAo, t — > 00. It is worth observing that the leading term in the r.h.s. of (|4.41|) 
does not depend on the regularization parameter e. It is also should be mentioned that equation 
(|4.41|) for the asymptotics of the Fredholm determinant exactly coincides with the one obtained 
in 1^ by independent and more direct methods. This fact partially justifies our, mostly formal, 
manipulations with the operator-valued Riemann-Hilbert problems. 



5 Positive chemical potential 

In the case of positive chemical potential the method of factorization, considered in the section 
^ is not applicable directly to the jump matrix ( 3.11| ). Indeed, now the determinants of the 
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operators G'ii(A) and G'22(A) have zeros on the real axis, and the scalar RHP (e) for the function 
A(A) = det ^(A) is ill posed. The factorization ( [4.22|) also does not exist, since the operators Q, 
P, P and Q have singularities on the real axis (see (|4.17|) - ( [4.21|) ). 

In this situation one can use the approach similar to the one proposed for the free fermionic 



limit of the problem in p9[. First, we have to specify the position of the determinants roots relative 
to the saddle point Aq. Denote the roots of the function, 

det Gn(A) = 1 - i9(A)Z(A, A)e<^^(^) = 1 - ^9(A)(1 + e-'*'^^^), 

as Aj, j = 1, 2, Ai < A2. We shall concentrate our attention on the case, 

Ai < A2 < Ao, 

assuming simultaniously that the roots of det 622 (A) lie to the left of Aq as well (see the arguments 
in the end of the section All the other possible cases of the roots position can be considered 



in a similar fashion. Moreover, similar to the free fermion case [^, the difference in the way how 
the roots are distributed along the real line does not affect the final answer for the asymptotics. 
What is important is the assumption that non of the roots is close to the saddle point. 

Secondly, observe that the integral operator I + V can be continued into some vicinity of the 
real axis, hence the integration contour can be slightly deformed in this vicinity. The Fredholm 
determinant evidently does not depend on such deformation. 

The deformed contour F is shown on the Fig.0. We will see below that this particular choice 




Figure 2: Deformation of the contour 



of the contour deformation will guarantee the solvability of the deformed scalar problem ( [4. 341) for 
the function A (A). 

As in the case of negative chemical potential we make the substitution 

X(A) = <I(A)7^(A), (5.1) 
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where as before 7?.(A) is diagonal matrix 

. ( ^(A) 

7^ = 



(5.2) 



V (^^(A))-^ 
and ^(A) satisfies the RHP 

1- ^. ^(A) ^ A — s> oo, 

2- ^. ^(A) is analytical function of A if A ^ F, (5.3) 

3^. ^_(A) = ^+(A) (^^(Ao - \)G^i + ^(A - Ao)(G'^2)"') , A G T, 

The difference between the Riemann-Hilbert problem (6), considered in the previous section, and 
RHP (/) is that the last one is formulated on the jump contour F instead of the real axis. Similary, 
both the operator-valued RHP for $(A) and the scalar RHP for A(A) = det ^(A) are set now on 
the contour P. Due to the choice made for the contour P, the function 

argdet (^(Ao - A)Gn(A) + 0(A - \^){Gl^)-\X)) 

can be defined in such way that it is continuous for all A 7^ Aq and approaches as A — ±00. In 
other words, the A-RHP has no index and hence is well posed. Using again the general formula 
(|2.14|) (with the real line replaced by P), we have that 



sign(5R^t-Ao) 



(5.4) 



(cf. (E33)). 

The jump matrix is given by the same equations (4.16)- ( [4.21 ) as before. The operators 
Q and P have singularities at A1.2, and hence the obstructions to the analytic continuation of 
the operator matrices M+ and M_ occure (the singularities of N± are not relevant since we have 
assumed that the zeros of det G22 lie to the left of Aq). Therefore, before factorizing the $- jump 
matrix let us make one more substitution (cf. ||29| ): 




$(A) = (A/ + fl)$(°)(A)(^^^ . (5.5) 

Here H is an operator satisfying the following conditions: 

(A2/ + n)<|(°HA2j 



(Ai/ + fl)<|(°)(Ai 




(5.6) 
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The RHP for the operator <I>'^'^)(A) has the form 

15. I>(°)(A) A^oo, 

2^. $(°)(A) is analytical function of A if A ^ i?, 

33. |)W(A) = <lf(A)G,(A), XeR, 

where the jump matrices Gg and (5$ ( |4.7| ) are related by 

' (A - Aa)? 



(5.7) 



(A - Aa)? 
(A - Ai) ■ i 







(A-Ai)z 



(5. 



Of course, now one should choose for G$ the operator g satisfying RHP (/). 

Since the roots Ai and A2 as well as the roots of detG'22(A) all are smaller then Aq, the 
factorization of the matrix Gg in the domain to the right from the saddle point produces the 
matrices N± whose analyticity in the relevant half planes is obvious. In the domain A < Aq the 
jump matrix Gg has the form 



Gs{X) 



Here 



'X-A^^ 
X-A, 



g«(A) = g(A) 



'A-A^' 
A-Ai 



(5.9) 



(5.10) 



and operators Q and P are given by ( ^.17| ) and ( ^.18| ). Thus, the jump matrix Gg can be factorized 
in the domain A < Aq as 



Gg = M+M_, 



where 







1 I P(i)e- \ 




M„ = 











(5.11) 



(5.12) 



Using explicit formulae ( ^.11 ) for the operators Gjk one can write 

g«(A)=g«(A) {f^)'\x)\2){l\rA>^\ 
pa)(A)=p«(A)^_(A)|l)(2|^^(A), 

^ 27rzZ(A,A)(^(A)-l)e-/'W ( X-A{ 
^ ^ ' I- Z{X,X)d{X)e't'oW \X-A2, 



where 



(5.13) 
(5.14) 

(5.15) 
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Note that the denominator in these formula is just the determinant of the operator Gu: 

det Gii(A) = 1 - Z{X, X)^{X)e'^^^^\ (5.17) 

and it has zeros when A = Ai 2. However, it is easy to see that p^^\X) has no pole at A = Ai as 
well as ^''■^■'(A) has no pole at A = A2. Thus the function p*^^)(A) (and, hence, the matrix M_) can 
be analytically continued into some vicinity of the real axis to the lower half-plane. Similarly the 
matrix M+ can be continued into a vicinity of the real axis to the upper half-plane. Hence, we 
again arrive at the contour Fig.[l|. 

The following considerations are absolutely identical to the considerations of the previous 
section. Introducing the operator U{X) by equations ( |4.25| ) ( with $ replaced by ^^^^), we again 
obtain that all the jump matrices associated with the function U{X) are exponentially small for 
the large value of t on the contour Fig.|l] away from Aq. Therefore, we have again that 

U^i. (5.18) 

However, in distinction of the previous case, in order to find the original solution x we need to 
take into account the contribution of the operator H ( ^3|) . 
The system of equations (|5.6| ) gives 



Ilii = -Aiz + Ai2<l>S^(A2)D, 

ni2 = A2il>S^(A2)D |.S?(Ai) (<|.g(Ai))"\ 

= Ai2l'^?(A2)D, 

1122 = -Aiz + A2l<l>S^(A2)^ <l>S°2'(Ai) (^^(Ai 



(5.19) 



-1 



where Ajk = A^ — Aj. and 



D 



dff (A2) - ^^^(Ai) (<l>g)(Ai))"' 1.5? (A2) 



(5.20) 



The solution <I>(''^(A) of the RHP [g) is asymptoticly equal to M_(_^ and M_ in the domains ^2 
and ^3 respectively. Hence, 

|.W(Ai)^M-(Ai)= I I ^^'^(^1)^^^^^^ 1, (5.21) 
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Thus, we obtain 



.(0) 



(0), 



I 

-g(i)(A2)e-^(^2) i 



<I>5°)(A0, <fi2^(A2)^0, 



(0), 



$(?(A2)^-g^^HA2)e 



^^"^(Ai) ^P«(Ai)e^(^i^ 
Therefore, equations ( |5.19| ) can be replaced by the equations, 

fin ^ -All + Ai2^, fli2 ^ A21L' $S(Ai), 



n21 ^ Ai2$g^(A2)/^, n22 ~ -Ai? + A2l<f^^AA2)/^ <fr2^(Ai) 



(0), 



.(0), 



and 



D 



?-<I>S°2^(Ai)|.(°HA2 



-1 



Define new vectors 



|1) 



5,{u- Ai)Z{u,Ai) 



|2) 



^(A,,A, 



-Z{u,Ki 



(1| 
(2 1 



^(A„A, 



\ M(A.) 

'AC(Ai)^(A,,,A,)' 



^ M(A.) 

where z = 1,2, and introduce two quantities 

a= (l|f?;i(A2)f?-(Ai)|i), 

/3= (2|f,^(Ai)(f.^(A2))"'|2). 

Then it is easy to see that 

$g^(A2) ^ -g«(A2)e-^(^^) (^^(A2))"' |2)(1| ^;^(A2), 

|.S?(A0^p(^)(A0e-(^^)^_(A0|l)(2|^^(Ai), 



and hence 



/3p(^)(Ai)g(i)(A2)e"(^i)-"(^2) 



1 + a 



/3pW(A,)g(^)(A2)e-(AO-^(A.) ^-(^^)'^)<^'^"(^^)' 
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n 



11 



-^2^ + 



n 



12 



1 + a/5p(i)(Ai)g(i)(A2)e-(^i)-(^2) 
Now, let us imtroduce the coefficients of the asymptotic expansion 

$(o)(A)=/ + l5W + 1(7(0)..., 
A A 



(5.31) 
(5.32) 

(5.33) 



Since ^^^^ ~ / for complex A, we have that B^^\C^^^ ~ 0. Using expansion ( [5. 331 ) and expansion 
(|4.9|) for the function g, we obtain the following asymptotic representations for the operators hkj 
and Ckji 

611 = fill + ^0 + As? + 0(1), (5.34) 

, 2 



C22 - cii = (Aif + n22)(Ai? - qI) - (Aaf + flii)(A2f + ^0) - Qi - fi + (qo) + ^(1)' (5-35) 



and 



612 =fll2 + 0(l). 



(5.36) 



The logarithmic derivatives of the Fredholm determinant with respect to x and t (i.e. Aq and t) 
are equal to the traces of the operators ( |5.34|) and (|5.35|) respectively. Consider more detailed the 
first one. From ( |5.31| ) and ( [5.34| ) we have 



&11 = ^0 + 



hence 



trfoii = tr^o + 



/3A2ip(^) (Ai)g(^) (A2)e"(^^)-"(^^) 
1 + a/3p(i)(Ai)g(i)(A2)e-(Ai)— 

a/3A2ip(^) (Ai)g(i) (A2)e"(^i)-"(^2) 



^_(Ai)|l)(l|^;i(A2) + o(l) 



(5.37) 



1 + a/5p(i)(Ai)g(i)(A2)e-(^i)— (^2) 



0(1) 



tr 



d 



go - ^^ In (l + a/3p(i)(Ai)g(i)(A2)e^(^^)-^(^^)) + o(l). 



(5.38) 



The trace of the operator ^0 as before is equal to the ffist coefficient of the asymptotic expansion 
of its determinant Aq (see ( |4.11|) ), which due to (|5.4| ) is equal to 



1 /" 

Ao = — dfi sign(Ao - ^jj) In (l - i9(/i) (l + e 



4>{fi) sign{R^-A(); 



(5.39) 



In comparison with ( [4.3^ ) the difference is that one should integrate now with respect to deformed 
contour F instead of the real axis. 

Similarly one can find the trace of the operator C22 — cn. We do not present here the details of 
these calculations; they are quite strightforward, although rather long. The resulting formula for 
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the leading term of the asymptotics of the Fredholm determinant is 
det{i+V) = (l +«/?p«(Ai)g(i)(A2)e"('^^^-"(^^)) 

X exp 1^ ^ dfi\x - 2/it| In (l - (l + ^^^{^i^-M)^^ . (5.40) 

The symbol \x — 2/i)f:| is understood according to the equation, 

\x — 2jit\ = [x — 2 fit) sign(3fJ/i — Aq). 

We would like to draw attention of the reader to the difference between the cases of positive 
and negative chemical potential. In the last case the determinants of the operators Gu and G22 

det Gii(22) (A) = 1 - ^(A) (1 + e- W'^(^)) (5.41) 

have no zeros at the real axis, and the zero-index condition, 

argdet(^(Ao-A)G'ii(A) + 0(A-Ao)(G'^2)"^(A)) ^0, |A| ^ 00, (5.42) 

takes place on the real line. Thus the integral in ( [4.41| ) is well defined. In the case of positive 
chemical potential det Gu have real roots, and therefore one should choose the integration contour 
more accurately in order to satisfy (|5.42| ). We have seen, that this means that the real axis have 



to be replaced by the contour T. 

The second important difference, is that the leading term of the solution of the RHP (a) is 
diagonal in the case of negative chemical potential ( [4.31| ). Thus the coefficient 612 is asymptoticly 
equal to zero (in the next section we will see that actually 612 = 0{t^^^'^)). This means that in 
the case of negative chemical potential the accuracy we reached so far is not enough for evaluation 
of the full combination / dudvbu ■ det (^1 + V^, which describes the correlation function. On the 
contrary in the case of positive chemical potential we obtain non-zero value for 612 already in the 
leading term of the asymptotics, 

- 1 + a/.p...(A.).' .;(i,),.,A.,-.,A„ ^^-(A0|l>(2| g(A.) (5.43) 

+ 0(1). 

Observe, that the pre-exponent factor in ( ^.40|) cancels the denominator in (|5.43| ), so we obtain 

/OO 11 
dudvg4A,)\l){2\f_iA,) 
-00 

X exp 1^ ^ dfi\x - 2fit\ In (l - ?9(/i)(l + e<^M«is>^{«M~^o))) | . (5.44) 

As we have mentioned already, the different asymptotics for negative and positive chemical 
potential has deep physical origin. This difference becomes essentially important in the low- 
temperature limit. 
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6 The localized RHP 



The estimates of the previous sections are not enough for comprehensive analysis of the correlation 
function ( p..l[ ). In particular, we saw that for the case of negative chemical potential the leading 
term of the asymptotics of the RHP (a) solution is diagonal, therefore we can not find the factor 
bi2, which is necessary for study of the correlation function. In order to improve the estimates 
obtained in the previous sections one needs to perform the second step of the nonlinear steepest 
descent method, i.e. to consider the RHP in the vicinity of the saddle point — so called localized 
Riemann-Hilbert problem. Setting f/ ^ / in ( [4.31| ) we neglect the contribution of the vicinity 
|A — Ao| < Now we consider this problem in more detail. The difference between positive 

and negative chemical potential now is not very essential, therefore for simplicity we concentrate 
our attention at the negative chemical potential only. 

Since the vicinity of the saddle point |A — Ao| < t~^^'^ is small for large t one could replace 
jump matrices N±{X) and M±(A) by their values in the point Aq: iV'±(Ao) and M±(Ao). However, 
N±{X) and M±{\) are well defined only in the vicinity of Aq, but not exactly in the saddle point. 
Indeed, all of them depend on g (see ( CT) , ( CT) , (CT) , (J^^, (lO^ ) and (WM) )- In turn, g 
satisfies RHP (b) with jump operator having a discontinuity in the saddle point. Recall the jump 
condition for the boundary values of g: 



where 



It is easy to see that 



g_ = g^V \e R. 
V = 9{Xo - X)Gn + ^^(A - Ao)(G^2)"'- 



hm P(A)^P,= hm I?(A). 

A — ^Aq— U A — ^Aq+U 



Thus, the solution g of the RHP (6) has a branch point at A = Aq. 

The explicit expression for jump operator P is given by (|3.11 
structure 

V = 9{Xo - A)(z + MXp) + 9{X - Ao)(z + ^(A).:;). 
Here /i(A) and /2(A) are equal to 

/i(A) = -^(A)Z(A,A)e^-W, 



(6.1) 

(6.2) 
(6.3) 

It has the following 
(6.4) 



/2(A) 



t^(A)Z(A,A)e'^-^W 
1- Z{X,X)^{X)eM>^y 



(6.5) 



and 



cI;(A) = Lu{X\u,v) = |1)(1|. 



(6.6) 
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The operator Cj will play an important role. Obviously Cj is one- dimensional projector with the 
trace equal to 1, hence, 

Cj"" = Cj. (6.7) 

Defining an operator-valued function (A — Aq)^'^, where ^ is a number, as a formal Taylor series 
with respect to Cj we obtain 

(A - Ao)«^ = i + uj({X- Xof - l) (6.8) 

We shall consider formula ( |6.8[ ) as non-formal definition of the operator- valued function (A — Aq)^'^. 
We will also assume that the branch of the function (A — Aq)^ is fixed by the condition, 

- TT < arg(A - Ao) < tt. (6.9) 

It is easy to check the following properties 

((A - Ao)«^)"' = (A - Ao)-«^ = t + Cj({\- Ao)-« - l) , (6.10) 

(A - Ao)f (A - Ao):^" = I + ^(Ao - \)Cj (e^^'^ - l) = 0(Ao - A)e2"^«^ + ^(A - Ao)f, (6.11) 
Put 

^(A) = ^('=)(A)(A-Ao)^^^°, (6.12) 
where Cjq = uj{\q) and s is a complex parameter. Then jump condition for g^'^^ takes the form 

Q^-^ = ^?^o, (6.13) 

where 

VoiX) = (A - Ao)V'^»^(A)(A - Ao):'^''°. (6.14) 

We want equation (|6.12|) to represent the singularity of ^(A) near the saddle point Aq, i.e. we want 
the ±-values of the factor ^'■'^•'(A) to be continuous at Aq. To ensure this property we demand that 
'DqW has no jump at Aq (cf. the analysis of the usual RHPs near the points of the discontinuities 
of the jump coefficients [0]). Observe that the limit values of i'(A) in the saddle point Vi and Vr 
commute with (A — Aq)*'^'^". Thus, using (|6.11|) we have 

lim Vo{\) = e-^^'''^'ii + MXo)C;o), 

(6-15) 

lim r'o(A) = i + /2(Ao)t2'o- 

A^Aq+O 

Making right-hand sides of ( |6.15| ) equal to each other we find 

1 l + /i(Ao) 

^TT 1 + /2(Aoj 
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Figure 3: Jump contour and jump matrices in the vicinity of the saddle point 



or, using (|6.5|) 



27r 



"(l - ^9o^oe'^°^^"^) (l - i9o^e'^^(^o))] , (6.17) 
where 

^9o = ^9(Ao), Zo = Z(Ao,Ao). (6.18) 
Thus, we have found the behavior of g in the saddle point. It is given by (|6.12| ), where g^'^^ has 



the limits as A approaches Aq from the upper and from the lower half planes, operator uj is defined 
in (|6.6D, s is given by ( |6.17| ). 

Let us turn back now to the RHP [d) (|4.26|) for the operator U. We consider this problem 
in the vicinity of Aq. The corresponding jump contour is shown on the Fig.^. It consists of four 
branches. One can continue the contour on the Fig.Q from the vicinity |A — Ao| < t~^^'^ to infinity 
and consider new RHP with jump conditions on two infinite axes. Such a replacement, of course, 
gives a non-zero contribution in the solution, but this contribution is exponentially small due to the 
presence of factors e^"^ in the jump matrices (cf. the classical analysis of the oscillatory integrals). 

Let us introduce a notation VV, which will denote matrices M_|_, M_, or for the corre- 
sponding branches of the contour. Due to the fact that the functions ^±^(A) are continuous at Aq, 
we replace W by Wo = {M±\ N±}, where 



igi'^^)-\X - Ao)-^^^" ■ g(°) ■ (A - Ao)-^^^«(^fV'-^ 



iV^^= I , I- (6.20) 



U'^^r\\ - Ao)-^^^" • • (A - Ao)-*^^«(^i'=V'e- 
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(0) 



I 



(6.21) 



(0) 



£•0 



(A - Ao)'""° ■ P(°) ■ (A - Ao)'^"'' 



(6.22) 



Here ^o"^"* = ^^'^■'(Aq), and 



g{o)g-r{Ao) ^ G2i(Ao)(Gii(Ao))-\ p{o)e-(^o) = (Gn(Ao))-^Gi2(Ao) 

(5(0)g-r(Ao) ^ (G22(Ao))-^G2l(Ao), P^e^^^") = Gl2(Ao)(G22(Ao))^^ 



(6.23) 
(6.24) 



Strictly speaking, one should add the subscript to ^q"^^ which appears in the matrices m'^\ 
N^^ and the subscript '_' to Qq"* which appears in the matrices M^^\ Unless the distinction 

is important, we usually will omit the subscripts to avoid the overcomplication of the notations. 

Based on the experience with the classical NLS (see |^ and [0), we expect that the replace- 
ment W Wo will enable us to capture the terms of order in the asymptotic solution of the 
original operator-valued RHP. These terms, as it has already been explained in the introduction, 
are of the most importance for the comprehensive asymptotic analysis of the correlation function 

O). 

Replacing W —>■ Wo, we approximate the exact solution U{X) by the operator-valued function 
U^^\X) which solves the following localized RHP (cf. again [^] and [^ ), 



00, 



l^^7(°HA)^/ A 

2''. U^°\X) is analytical function of A if A ^ C, 



(6.25) 



where C is the contour depicted in Fig.^ In this and the next two sections we will construct an 
explicit solution of this problem. 



Making the substitution 



where 



S 



(A-Ao)-'^'"°(^fVe-^(^)/2 Q \ 

(A-Ao)^^^''4^VW/^ J 



(6.26) 



(6.27) 
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we obtain new RHP for the operator ^: 

2\ ^(A) is analytical function of A if A ^ C U i?, 



where 



and 



M 



(0) 



M 



(0) 



£ 
Qo i 

t 



' +,'1' 



z 

? p(o) 

i 



(6.28) 
(6.29) 

(6.30) 

(6.31) 

(6.32) 



Since the matrix S depends on (A — Ao)*^'^° and ^0*^^ we have an additional jumps on the real axis 



3*2. ^_ = ^,K, XeR 



where 



/ 



K = K, 



K = K_ 



V 



? + /i(Ao)c2;o ^ 

(i + f^(X^)u;^)-^^ 

2 + /2(Ao)cI'o ^ 

{i + U\o)u;^r' ; 
For the functions /i,2(A) see (|6.5|). 



A < Ao, 



A > Ao, 



(6.33) 



(6.34) 



(6.35) 



It is convenient to make the following change of the variable A, 



X ^ z : A — An = 2;e*4 , 



so that, 



S 



V 







(6.36) 



Note that matrices VV,j, and K do not depend on A, hence they do not change. 
Let k be the first coefficient of the asymptotic expansion of U^^^: 

u^'^ = i + ^ + .... 

X 



(6.37) 
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Then the asymptotics of ^ is given by the expression 



K 



(6.38) 



Due to the property that all the jump matrices do not depend on z we conclude that (5^^)^"^ is 
holomorphic in the finite complex plane. The asymptotics of this 'logarithmic derivative' (S^^)^"^ 
can be found from equality 



It is easy to see that 



/ ^ is 
I tzi 



V 







IS 



-tzi H cJq J 



therefore 



Due to the Liouville theorem (S^^)^^^ is a first order polynomial. Using ( |6.38| ) we find 

= -tza3-t[k,cr3]e-'^, 
and, hence, ^ satisfies ordinary differential equation 

a,^ = -t(za3 + [/i,cr3]e-'?)^. 



(6.39) 



(6.40) 



(6.41) 



(6.42) 



(6.43) 



The arguments which led us to the differential equation ( |6.43|) are identical to the arguments 



used in the theory of classical integrable systems. Similar to the pure matrix case p9[, [B^], El 



the first order 'matrix' differential equation for \l/ implies the second order differential equation for 
components ^jfc. After one more replacement, ^ = zy/2i, we obtain, for example, for 



11 



2 4 



11 



z>^ii = 0. 



(6.44) 



Here v = —2itki2k2i. This equation looks like the parabolic cylinder equation, except that v is 
an operator. In the next section we shall give rigorous definition of solution of equation ( |6.44|) . At 
the moment one can consider the parabolic cylinder function (PCF) D(,{C,) with operator- valued 
index z> as a formal Tailor series with respect to 0. This series satisfies recurrence, which is valid 
for usual PCF: 

2 



(6.45) 



0, 
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Using these properties, one can check directly that the following operator-valued matrix 

^ = I \L. (6.46) 



satisfies equation ( |6.43| ). Here 

P = — 2itKi2K2i, z> = 2itk2iki2, (6.47) 
and L is piecewise constant operator-valued matrix. 

7 The parabolic cylinder functions with an operator index 



To complete the solution of the model RHP (h), or equivalently the RHP (i) (see ( |6.28| ) - ( |6.32| )), 
the quantaties kji and L involved in (|6.46| ) are needed to be determined explicitly in terms of the 



jump matrix VV*. Also, it follows from the comparision of expansion ( |6.37| ) for and expansion 
(|2.16|) for original operator x that 

ki2 = bi2 and /t2i = &2i- (7.1) 

Thus, in order to obtain the asymptotics of the correlation function we need to find just these 
operators. 

To achieve the objectives formulated above, we have to make the PCF Do{C,) with the operator- 
valued index a non-formal object. In fact, we need to control its asymptotic behavior as ^ ^ oo 
everywhere on the complex plane ^. In other words, similar to the classical PCF, we want -Dj>(0 
to be a genuine special function. To this end we need to know more about some special features of 
the operators kji and z> which follow from the setting of the RHP(z). In this and the next sections 
we show that the structure of the jump matrices of the model RHP(i) suggests a very special 
anzats for these operators which allows eventually to reduce the operator-valued problem {i) to a 
matrix one. This in turn will give us the possibility to obtain the explicit expressions for kji. 

In the previous section we have already replaced variable A by variable ^: 

A-Ao^-^e-f. (7.2) 



Further we shall work with variable ^ only. Therefore we present here jump contours in the complex 
plane of ^. The complex ^-plane consists of six sectors /,..., VI. The boundaries of the sectors 
Ci-j, i,j = VI are shown on the Fig.^. They are jump contours of the matrix The 
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Figure 4: Jump contour and jump matrices in the -plane 

corresponding operator- valued jump matrices also placed on the Fig.^. In the parenthesis we have 
placed scalar analogs of jump matrices, which will be discussed later. 

The piecewise constant matrix L and operators ki2 and K21 are uniquely defined by the condition 
that matrix operator ^ ( |6.46| ) satisfies the jump equations indicated in ( |6.28| )-( [0^ ) and possesses 



asymptotics (|6.38|) . These restrictions should provide us the explicit expressions for kjk and L. 

In what follows we shall consider vectors |1), |2), (1| and (2| only in the point A = Aq. Therefore, 
we introduce new objects 



° (5, M- Aq Z n,Ao o Zq 

K{\o)Zo V M Ao 



|2) = Jt7?^^(-,Ao), (2| ^^(--^o)^(-,Ao) 



(7.3) 

Recall once more that Zq = Z{\q, Aq). Denoting as before (u)(Ao) = tJo, we have 

u;o=\l){l\, u;^=\2){2\. (7.4) 

Using the Taylor series definition, we obtain for an arbitrary holomorphic at the origin function 
fix): 

=/(o) (£-|i)(i|) + |i)(i|/(x), 

^ ^ ^ ^ (7.5) 

/(x|2){2|) =/(0) (?-|2)(2|) + |2)(2|/(a;). 



00 / o o \ 00 

|2)(l|/(x|l)(l|)=/(a:)|2)(l| 
|i)(2|/(x|2)(2|)=/(a:)|i)(2| 



(7.6) 
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Note also the following generalization of ( [7. 51 ), 

/ (ai + ) = /(a) (z - ) + fix + a), 

f {at + x|2)(2| ) = /(a) (z - |2)(2| ) + |2)(2| f{x + a), 

which takes place for a function f{x) holomorphic at the point x = a. 

Now, we are ready to describe our anzats. The structure of jump matrices N±^q, and 
allows us to look for the operators ki2 and in the form of one-dimensional projectors. Therefore 
we set 

ku = 42Qii\hmQl^r, (7.8) 
«:2i = 4i(^if)"'|2)(i|(^a)-^ (7.9) 

Here fi:^2 nfi are constants, which should be defined. One should choose in both formula the 
sign plus for sectors J, II, VI, but the sign minus for sectors ///, IV, V. However, operators kjk, 
being the coefficients of the asymptotic expansion of the operator U, must be the same in any 
sector of complex ^-plane (since, by definition, they do not depend on ^). Therefore, the values of 
K^j^ and K~j^ on the contour C//_/// U Cv~vi should be related with the jump of the operator Qq"^ 
on the same contour, in order to provide the unique value of the operators kjk. Thus, we demand 

KtJoimi gif = KiJ,^l\l){2\ gi^l^. (7.10) 



Using (|03|) - (|6J[5|) we find 



Uiy'il = {i+f2iXo)\my' 



(7.11) 



This leads us to the jump condition for k,i2 

'^12 

Similar consideration gives 



42(l + /2(Ao))-'. (7.12) 



'^2'i = 4i(l + /2(Ao))'. (7.13) 
Substituting (|7.8|) , (|7.9|) into expression (|6.47|) for z> we obtain 

u = ugi^^\l){l\igY\ (7.14) 

and 

^ = -K^r)-'|2)(2|4^^^ = -z>^, (7.15) 

where 

l> = —2itK,i2K,21 — —2ittii2^2i- (7-16) 
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It is worth mentioned that due to ( |7.12| ) and (7.13) the number v (not only the operator z/ !) is 
the same in all sectors of ^-plane. 

Using explicit expressions for the operators v and v ([7. 14| ) , ( [7.15[ ) and formute ( |7.5| ), ( [7.7| ) one 
can define 'operator-indexed' PCFs as 



D,{i) = DoiO [i - |1)(1| ) + |1)(1| D^iO 



(7.17) 



D,-dO = Qo^ D^iiO - ) + |1)(1| /^.-i(0 



'AY'. 



DM) = iin'' [Dom - |2)(2| ) + |2)(2| D_M) 



4c)T 

Qo ■ 



(7.19) 



D^-,m = AY^' [D-M) {i- |2)(2|) + \2)(2\D^^^^{^0 



4c)T 

Qo ■ 



(7.20) 



We would like to mention that u in these formulas is a number, and therefore Dy{S,) is the usual 
(not operator- indexed!) PCF. It follows from ( |7.11| ) that 



(7.21) 



and 



(4^i-)-|2)(2|^r = (^r)-^|2)(2|^if. 



(7.22) 



Therefore, the r.h.s.s of ( [7.17| )- ( [7.20| ) do not have jumps on the boundaries Cn~iii and Cy-vi- 
One can check once more, that functions, defined by (|7.17| )- (|7.20|) satisfy recurrence (|6.45|) . 



8 Reduction to the matrix RHP 



Define two mappings Aq{M) ('parallel mapping') and A^{M) ('orthogonal mapping') of 2 x 2 
matrices into algebra of 'hat '-operators. The first mapping is defined for arbitrary matrix M: 



iS(M) 



^ Mn|i)(l| Mi2ri)(2| 
^ M2i|2)(i| M22|2)(2| 



The second mapping is defined only for diagonal 2x2 matrices: 

/ (?-|i)(i|)Mn 



|2)(2|)M22 



12) 
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Similar mappings were introduced by Korepin, and they have aheady been used in [^. Similar 
to 1^ it is strightforward to check that 

M^Ao{M)=Ai{I)+Al{M), (8.3) 

where I is the unite matrix, is a representation of Gl{2, C) (cf. (|1.19|) , (|1.21|) in the introduction) 
. We would like to emphasize that in distinction of the representation, considered in [^, the 
mappings v4.o(M) and Aq{M) preserve the analytical structure of a matrix function M(A) since 

o o o o 

the projectors |1)(1| , |1)(2| etc. do not depend on A. 

The mappings ^o(M) and Aq{M) are not representations. In particular ^o(/) 7^ / and 
^0 (-^) 7^ ^- Nevertheless these mappings possess very important properties: 

Jl{M)Ai{N)=Ai{N)J^i{M) 



0, 



4(M)i|l(iV)=ill(MiV), 
Ai{M)Ai{N)=Ai{MN). 



•4) 



These properties allows one to reduce the operator-valued RHP {i) for ^ to two matrix RHP for 
parallel and orthogonal parts of its pre-image. 
Let 



7^(0) 



^ 9^' 



(c)Tx-l 



V (^D 



Since Qq has two values ^0+ and ^o_, the operator- valued matrix T^*^"-* also has two values:7^+ in 
sectors J, JJ, VI, and T^^"* in sectors JJJ, IV, V . The jump condition is 



n 



(0) 



where K_ is defined in ( |6.35| ). 

Let us introduce also two scalar matrices and 

\ Do{ii) 

( 



17) 



The matrix /3" depends on /ti2 and K21, which have two values in different sectors of the complex 
plane. Hence, /5" also has two values: /3| in sectors /, //, VI, and /3i in sectors ///, IV, V. The 
jump condition is (cf. (|7.12| ), (|7.13| )) 

/3i = (fc_)-Wfc_, (8.9) 
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where 

l + /2(Ao) \ 

(i + /2(Ao))-' y 
Let matrix-operator L has the structure 

where 

Ui, 



and 



V 







'22 



'11 ti2 



V «21 «22 



Then combining Eqs. (|7T7D -(f7:20 D , 
obtain 

It is easy to check, that the matrix 5*^^ (see (|6.36|) ) also has the structure: 
where 



/ 



f2 it>.i 



e 4 



p2 it\^, 



( 



e 4 



V 



e 4 



' 2 e~ 4 ^ 



V 



(8.10) 



fS.lll 



(8.12) 



(8.13) 



and using the properties of mappings ( ^.4|) we 

(8.14) 



(8.15) 



(8.16) 



(8.17) 



(2t)~e-+ 2 j 

Here we have used formula ( |7.4| ) for Cjq and o;^. 

Consider now the mappings for jump matrices. The jump conditions for the operator \1/ at 
the contours Cj-vi, Ci-.ii, Cni-iv and Cjv-v are given by the matrix VV^ = {iVj°^^, M±\} (see 
Eqs.(|]2|)-(|]32D and Fig.|). Formulae for Q^^\ P^^\ Q(°) and p(°) are given in (|]2|)-(|]2|). 
Using explicit expressions for Gjk (see ( p.ll|) ) we find 



Q(°) = gW|2)(l|, Q(°)=g(°)|2)(l|, 

p(°)=p(°)|i)(2|, p(°)=p(°)|i)(2|, 



^.18) 
^.19) 
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where 



,{0) 



p 



(0) 



1 - Zot9oe'^°(^«) ' 



Obviously, all the jump matrices VV<i- can be written in the form 



Here 



and 



m = AiiI)+Al{w^). 



= {n±,m±}, 



/ 1 
\^ 1 
/ 1 



n_ 



1 



1 

qW 1 

1 

1 



(8.20) 
(8.21) 

(8.22) 
(8.23) 
(8.24) 
(8.25) 



The jump condition at the contours Cji^ui and Cy-vi are given by matrix K±, which also can 
be presented as 



where 



(8.26) 
(8.27) 



(8.28) 



l + /i(Ao) ^ 

(i + /i(Ao))-^ ^ 
( l + /2(Ao) 

(i + /2(Ao))-^ 

The jump matrices m±, n± and k± are placed on the Fig.| in the parenthesis. 

Thus all operators entering the RHP (i) and its solution are written as the ^o()- mappings of 
2x2 matrices. Let us summarize the preliminary results. 

Using the mappings v4.q (•) and .Ao(-) we have presented the operator- valued matrix ^ as a sum 
of orthogonal and parallel parts (see ( ^.14j )): 



V 



^ = 7^w (^i^(/3^£^) + 4(/?"£")y 

We have similar representation for the operator 5*"^ 

5-l = 7^(°) (i^(s^) + 4(^" 



(8.29) 



(8.30) 
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This gives us the asymptotic conditions for orthogonal and parallel parts of ^. Indeed, since 
^ S~^, when ^ goes to infinity, then the pre- images of the parallel and orthogonal parts of 
should go to the pre-images of the parallel and orthogonal parts of respectively: 



Finally, all the jump matrices also are written in terms of orthogonal and parallel parts: 



(8.31) 
(8.32) 



=A^{I)+AI{w^), (8.33) 

k^ = A^{I)+Al{k±). (8.34) 

The parallel and orthogonal parts of ^ should separately satisfy jump conditions. Thus, we have 
reduced the operator-valued RHP to two matrix ones. Namely, we need now to find piecewise 
constant matrices £" and constants and such that: 

1. both of the matrices p-^i-^ and possess the asymptotics prescribed; 

2. both of the matrices jJ-^i-^ and possess the jump conditions prescribed. 



We start with the asymptotic condition for the orthogonal part. One should find matrix 
such that 

5— >oo 



e 4 



2 



V 







e 4 





2 it\ 



8.35) 



Here p-^ is given by (p.7|): 



P' 




(8.36) 



Since -Do(0 = 6xp(— ^^/4), we conclude that 



itXl 
expj ^0-3). 



(8.37) 



Let us discuss now the jump conditions for orthogonal part. Observe that orthogonal parts of 
all jump matrices are equal to Aq{I). Thus, matrix P^i-^ has no jumps on the contours C/.y/, 
Cj-ii, Cjii^iv and Cjv-v- It also has no jumps at the contours Cn-ni and Cy-vi in spite of 
orthogonal part of ^ depends on matrix iz^^\ which may have jump just at the contours Cn-iii 
and Cv-vi- Let us check the corresponding jump condition. We have (see 



-(c) 
^0- 



|1)(1| +(l + /2(Ao))|l)(l| 



f8.38) 
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(see ( |6.15|) ). Multiplying this equality by 2 — |1)(1| from the right, we obtain 



0^l(i-m\ 



(8.39) 



Thus, the orthogonal part of Qq^^ (and, hence, orthogonal part of '^^°'') has no jump at the contours 
Cn-iii and Cy-vi- 

The orthogonal part of \1/ is found: 



^.40) 



where and are given by ( |8.36|) and ( ^.37| ). 

Consider now the parallel part of ^. In order to simplify some formulae, we make replacement 



^exp(|-^ln(2<)-^-^^ 



(^3 



(8.41) 



First, one should find piecewise constant matrix £ such that matrix f3H possesses the following 
asymptotics 



Recall that /?" is equal to 



■ e 



(8.42) 



(8.43) 



The asymptotics of PCF D,y{^) strongly depends on argument of ^. In general form it may be 
written as 

^.(0 Ce-'^ + c{u)r'-'e^, (8.44) 

where constant c(z/) does not depends on |^|, but only on arg^ It easy to see that the 

asymptotics of (|8.42| ) can coincide with (|8.43|) only if 



ly = IS. 



^.45) 



Then we need to choose matrix C. in order to provide the correct asymptotics behavior in every 
sectors of the complex ^-plane. We drop out the detail of these calculations and present here the 
matrix £ in each sector: 



1 

Vo 



TT e 4 



t K2lT{u) 

e 2 J 



in the sector /, 



(8.46) 
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TT e 2 4 



\ \ t Ki2r(-Z/) 



TT e 



' 4 \ 



e 2 



in the sector //, 



(8.47) 



TT e 2 4 



V V t Ki2T{-u) 



t K2ir(z/) 



e 2 



in the sector ///, 



f8.48) 



/ 1 



TT e 2 4 



V V if: Ki2T{-u) 



\ 



e 2 



in the sector IV, 



.49) 



/ 1 



e 2 



in the sectors V and yj. 



(8.50) 



Let us illustrate how one can obtain these results. Consider sectors V and VI where — f < 
arg^ < 0. For these values of argument the asymptotics of the PCF is extremely simple: 



(8.51) 



D^M) e-- 



(8.52) 



Putting ii2 = ^21 = 0, in = 1 and £22 = we satisfy the asymptotic condition. It is worth 
mentioned that matrix i has no jump at the contour Cy-vi^ which is in completely agreement 
with the jump condition ( ^.651) below. Thus, we arrive at ( ^.50|) . 

Consider now sector I. Here the asymptotics of -0^(0 is just the same, but D^^ii^) behaves 
as follows 



Tiv) 



-e 2 



^.53) 



Putting again £21 = and in = 1, we obtain, that the first column of possesses correct 
asymptotics. But the second column has an additional term in comparison with sectors V and 

^ ^ _ ITT 

VI. In order to compensate this term, one should put £22 = and £12 = a/t-^i^- We obtain 
formula (|06|) . 

Similar considerations allow to find all other formula ( ^.47] )- (|8.49| ) for matrix £. 
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Finally, we should check, that all the jump conditions are valid. After replacement i the 
jump conditions for (3H have changed. Let us formulate them explicitly. 

The diagonal matrices kj- do not change under replacement ( p.41| ). However, the matrices m± 
and n± should be replaced by and respectively: 



m 



( ^ ° 
1 



1 p 



n_ 



1 




where 



P 



itXl 



(8.54) 

(8.55) 

(8.56) 
(8.57) 



hv = ^iiifn'_, 



Thus, the matrix £ should possess the following jumps: 

£vi = £in'_^, at the contour Cj-vi, (8.58) 

iim'_^_, at the contour Ci-n, (8.59) 

at the contour Cjn^iv, (8.60) 

iiv = ^y^'_, at the contour Cjv-v- (8.61) 

As for contours Cn^m and Cy-vii the situation is slightly more complicated. Recall jump 
conditions for iV-^^ and /?" 

7^L°^=7^f (ij[(/)+4(A:_)), (8.62) 

/5i = (A;_)~V1 ■ k_. (8.63) 
Then for the parallel part of ^ we have at the contour Cy-vi'- 

7tL°Uii(/?i£V) = 7^^U5(/9|fc_£V(A;-)-')4(^-) 



Hence, we find that 
Similarly 



= nfA{P\ivi)A{k-). 

ivi = k_iv{k-)~^ , at the contour Cy_y/. 
£jj = k_ijjj{k^)~^ , at the contour Cu-iii- 



(8.64) 

(8.65) 
(8.66) 
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One can easily check that the jumps on the contours Cy-vi and Cu^ni are vahd automatically. 
The jumps on the other contours give us explicit expressions for /tfg and nfi. Indeed, for example, 
due to ( ^.581 ) we have 



which implies 
and, hence. 



{il)l2 = -P 



n'(0) 



TT e 4 



^21 



tp'(o)r(z/)' 



Then due to ( |7.16| ) 



■'12 



2tK. 



(8.67) 
(8.68) 

(8.69) 
(8.70) 
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Using equalities ( |7.12| ) and ( [7.13| ) we can find k,i2 and k,2i- 

The jumps on the other contours give us just the same expressions for k^^. One should only 
take into account the identity 



].71) 



This identity can be checked directly, using explicit expressions for all variables entering (|8.71|) 
(see (p6D-(|33), 

Combining all necessary formulae after simple algebra we arrive at 



Kl2 



Here 



and 



7 = 27rZo(t9o - l)r(z/)(2t)-'^e^(^")+^+'^-**"§, 
1 



2m 



In 



0a(Ao) 



(8.72) 
(8.73) 

(8.74) 
(8.75) 



This gives us the solution ^ of the localized RHP (i). Let us present here once more the main 
formulae. The operator ^ is equal to 





(c)Tx-l 



V {gD 



(8.76) 
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The definition of the mappings (■) and Aq{-) is given in (^.11), ( ^.21) . The matrix is equal 
to 

' (e 



P^i^ = exp 
The parallel part of the matrix f3 is 




i.77) 



( 



/€2iy2te-^fD,_i(0 



(8.78) 



Here -Djy(0 is the PCF, v is given by (|8.75|) . The constants k\2 and K2\ have two values: for 
the sectors /, // and Vl^ and \^~^ for the sectors ///, IV and V (see Fig.|^). The values k'^^ are 
defined in ( |8.69| ), ( |8.7CI| ). The constants kJ^ are related with the last ones by ( |7.12D , ( |7.13| ). 



The matrix it is given by Eqs. (|^), (CT)-(|8:50|). 

For our purposes it is most important to know the operators kyi and k^x since the equations, 



kv2 = bi2, 

These operators are given by (|8.72|) , (|8.73| ). 



J21- 



(8.79) 



The case of the positive chemical potential can be delt with in a similar fashion. Instead of the 
operators Q*^"^ and P'-"-* one should consider operators 



= Q 



(0) 



'Aq- As' 
Ao-Ai. 



pm ^ p(0) 



An - A, 



(8.80) 



The corresponding replacement should be done for the scalars g^^-* and p^^^ as well. The rest of 
calculations is just the same as for the case of the negative chemical potential. The difference, 
however, is that for the negative chemical potential formulae (|8.79|) provide us with the leading 
term of the large time asymptotic expansion of the operators 612 and 621? while for the positive 
chemical potential the operators ki2 and ^21 define only corrections to the coefficients 612 and 
621- Indeed, in the last case, as we have seen, the leading term of the asymptotic expansion of 
bi2 is given by ( |5.43| ). In comparison with this expression, the operators ( ^.72| ), ( ^.73| ) have the 



order of hence, they are only corrections to the leading terms. Another words in the case of 
positive chemical potential the operators ki2 and K21 are equal to the coefficients of the asymptotic 
expansion ( ^.33|) b[^2 ^^cl B^i respectively: 



1^12 



^12 } 



1^21 



-D21 ■ 



^11 



Finally, we need to say few words about the operator ^q'^^ 
operator is defined by the equation. 



-(c) 
^0 



entering Eqs. ( PT2|) , (|T|). This 

(8.82) 
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where the operator- valued function ^'^'^•'(A) in turn is defined as a unique solution of the RHP 
(|6.13|) - (|6.14|) . This RHP is equivalent to the singular integral equation (cf. ( |1.15|) ), 



fi — X — iO 



(8.83) 



Therefore, to determine gf' one has to solve this integral equation. In other words, the time- 
independent operator factor ^g^'' in our final formulae is still a transcendental object. It is clear 
however, that it does not depend on the function V'(A) since the jump matrix T>q[X) only depends 
on the functions 0a(A), 0d(A), and the ratio Aq = x/2t (and, of course, on the chemical potential, 
coupling constant and temperature). We shall see later that this property is sufficient in order to 
calculate the asymptotics of the correlation function up to the numerical constant (depending on 
Aq, /i, c, and T). 



9 Differential equations 

In this section we continue to consider the case of negative chemical potential. 

The results of the previous sections provides us the leading terms of the large time asymptotics 
for the coefficients hjk- In particular, we have found that 



n2 



7Z/ 



(2t)-'^-i/2g^(Ao)-.ag^{c)|°i^^2|^{c)T^ 



J2l 



7 



(2t)''-i/V^("")+**^o(^ff )-i|2)(l|(4^i)-\ 



(9.1) 



(9.2) 



We have extracted explicitly the dependence on t and function ip{X). The constant 7 is equal to 

7 = 27rZo(t9o-l)r(z/)e— +-, (9.3) 
and it depends only on functions 0d(A) and 0a(A). Recall also, that 



In 



1 - ^9o^oe'^°^^"^) (1 - 79o^oe'^^^^°^ 



(9.4) 



Corrections to the asymptotic formulae (|0|), ( p.2| ) decay faster than t'^l'^ . However, in the poduct 
^12^21 the oscillations cancel, and after integrating equation ( |2.2(JD twice these corrections might 
produce the non- vanishing contributions into det(/ + V^). 

In order to make our asymptotic estimates more accurate, we use differential equations for hyi 



and 621 (Hi 



idth2i + 5^fe2i = 2621&12&21, 
- idthvi + f^x^i2 = 26i2&2i&i2- 



(9.5) 
(9.6) 
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We shall look for the solution of this system in the following form 

6i2 = (2t)-'^-^/V^*"oC(Ao,t), 



J21 



u-l/2 itXl 



57(Ao,t). 



(9.7) 
(9.8) 



It is convenient to rewrite the differential equations (|9.5| ), (|9.6| ) in terms of variables Aq and t. 
After simple algebra we arrive at the differential equations for the operators C, r): 



- -{CfiC - 

dt V ' ^ 
- jivCfj-iufj) = 



1 






{[ 













{u'Y ln\2t) - v" ln(2t) C - 2z/'C' ln(2t) + C" , 



'\2 1 2/ 



(9.9) 



(9.10) 



Here prime means the derivative with respect to Aq. It is easy to see that we can assume the 
asymptotic expansions pO] obtained for the solutions of the scalar Nonlinear Schrodinger equation: 



oo In 



C = Co(Ao) + EECn.(Ao) 



(ln(2t))^ 



n=l fc=0 
oo 2n 



^ = 57o(Ao) + E E'7nfc(Ao 



(ln(2t))^ 



(9.11) 
(9.12) 



n=l fc=0 



where coefficients and r)„fc depend on Aq only. Substituting these decompositions into (|9.9|) and 
(|9.10|) we obtain for the leading terms Co and r^o 



wC.Q = Co^oCo, 
ii^Vo = 570(0^70, 

These relations ( p.l3| ) are valid automatically, since due to Eqs. 



(9.13) 



^) we have 



Cofjo = iug''Ql\l){l\ {g'^Q 



In turn for the first corrections Cifc and r]ik we obtain the following recurrence 

Ci257o + C01712 = 0, 
Cu^o + Co^ii = ^ ('^'Co^o) , 
Cio^?o + Co^?io = f^'Co^o) + 7 (Co^7o - Co^o)' • 



(9.14) 
(9.15) 



(9.16) 



2i 



4 



Similar to the scalar Nonlinear Schrodinger equation the relations ( |9.16| ) allow us to find more 
precise estimate for the product = —idxbu. Recall that in order to obtain asymptotic 
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expression of the Fredholm determinant we need to find only trace of bu, but not bn it self. 
Obviously 

trCo'7o = ^i^- (9.17) 
With the trace of the expression (Co^o ~ Co^o) — the situation is more complicated, since we do 

^(c) 

not know explicit expression for . On the other hand it is possible to extract the dependence 
on the function '?/'(A). Using formulae ( |9.1| ) and (P^), we find 

CoV'o - Co^o = -2^z>^'(Ao) + Co, (9.18) 
where Cq does not depend on the function '?/'(A). Thus we have 

trfo'n = -u + ^(z/V)'(ln(2t) + 1) - ^(z/^'(Ao))' + ^ trCp + 0{\n't/t'). (9.19) 

Note. One could expect that due to the ( |9.11| ), ( |9.12D the last equation is valid up to the terms 
of order In^t/t^. These corrections do enter the expressions for b^ and 621 ? however they vanish 
(as well as the terms In^t/t^) for the combination tr 612621- 
Finally one should use the identity 

^aAolndet(/ + y) =ztr6n, (9.20) 



and after integrating of ( p.l9| ) with respect to Aq we obtain improved formula for Fredholm deter- 



minant. Of course, the integration constant may depend on t, therefore we need to use also the 
relations 

(^dt - y^Ao^ lndet(J + V) = ar (622 - Cn) , (9.21) 

and 

dtbu = ^ (5ao6i2 ■ hi - 612 ■ dxM) . (9.22) 

These equations allow us to find integrating constant up to some function, which might depend 
on Ao since we did not specify Cq in (|9.18| ). We skip the details; they are very similar to the 



free fermionic case spelled out in [^. Thus, we have for the Fredholm determinant the following 
estimate, 

lndet(J + r) = ^ r (\x-2Xt\+t sign(Ao - A)^^) 
xln{l - (l + e^(^) «is>^(^-^o)^ I 

-y ln(2t) +Co + C(ln2t/t). (9.23) 
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Here Cq includes integrating constant and other terms depending on functions 0^, 0d, the ratio 
Aq and the non-dynamical physical parameters h, c, and T. 

In order to calculate the correlation function of local fields we need to find the combination 

B = r dudvhi2{u,v) ■ deiii + V). (9.24) 

The formula we obtained for the Fredholm determinant and the operator 612 allow us to find the 
asymptotic expression for this object 



B = C(0^,0A)(2t)- 2 e'^(^«)-^*^o 



X 



ln{l - ^9(A) (1 + e*W-s'^(^-^o)^ | (i + q^i^2 ^^-^^ _ ^g ^g) 



Here we have extracted all dependence on t and function ipi^X). The factor C{(f>D, (Pa) depends no 
Aq, h, c, and T, and it does not depend on ^ip^X); parameter u is equal to 



1 



In 



(1 - i9o^e'^°(^°)) (1 - ^9o^oe"^^^^°^)] . (9.26) 



2iri 

and the function 0(A) is defined by the equation, 

0(A) = 0a(A)-0d(A). (9.27) 

In the case of positive chemical potential the calculations are quite similar. It is easy to show 
that the operators ^^'^ ( |5.33D satisfy the same system of the differential equations ( |9.5D , 



(|9.6|) as the operators 612 and 621- The asymptotic behavior of By^ and ^21^ was found in the 
previous section — Eq. ( p.81| ). Thus one can use the same approach in order to obtain more precise 



estimates. We do not present here the details since in the section 11 we will develop a modified 
method allowing to obtain the most important, i.e. ip'{\) - correction to the leading term of the 
asymptotics directly, without using the differential equations. 

Finally, we need to say several words about the regularization introduced in the section ^. The 
only object in ( |9.25|) , which depends on e, is the factor C{(f)D,4>A), since it depends on Qq"^ and 

o o 

|1)(2| . The exponential and power law factors do not depend on the regularization. We shall see 
later that C{(j)£,, 0a) is only responsible for a numerical factor in the asymptotics of the correlation 
function (Ol) 



55 



10 Modified approach. Preliminary discussions 



Assuming the dual fields tp{X), 0a(A) and to be classical functions, we found the asymptotic 

expression for B, which is the combination of the Fredholm determinant det (^1 + and factor bi2 
(|2.18|) . This asymptotic expression has the form 

where J^c, and J-'e are functionals of the dual fields. Their representations are given in ( |9.25|) . 
The functionals J-'e((j)A, 4>d) and J?p(^0A, ^d), defining the exponential and power laws respectively, 
are found explicitly 

iJ^e(0A, 0d) = -itXl + ^ 1^ - 2At| ln{l - {l + e<^W-g°(^-^o)) } d\ (10.2) 



^p(</)a,</>d) =-^(^^ + 1)', (10.3) 



where v is given by ( |9.26D . The functional J^c depends only on the ratio Aq = x/2t, which is fixed. 



hence it behaves as a constant when t goes to infinity. This functional consists of two parts, which 
play essentially different roles in the process of averaging with respect to auxiliary vacuum. Let 
us extract from JF^ the part, containing the field if), namely 

= ^cA4>A, 4>D, 'ip)^c4>{4>A, 4>d), (10.4) 

where 

J'c^{<Pa,<Pd,^) =^{Xo) + ^ J^jign{Xo-X)^^ (10.5) 

For the functional J-'ccf,, which does not depend on tp, we do not have any explicit formula; this 
functional depends on the operator go'^ Thus, the asymptotic behavior of B is found up to the 
constant factor JF^^. The similar factor in the free fermionic case can be determined by the use 
of the auxilary differential equations with respect to the temperature and chemical potential (cf. 
[pup. Here we do not pursue such goals, and we restrict ourselves by establishing the exponential 
and power asymptotic laws. Thus, in the case of the dual fields being classical functions, we have 
complited the analysis. It should be also mentioned that the exponential law which we found 
coincide with the one obtained in via the independed and more direct methods. 

However, as we have already mentioned in the introduction and section |], the operator nature 
of the dual fields strongly changes the situation at the stage of calculation of the vacuum mean 
value (0|i3|0). One must be sure that the corrections to the leading term do not contribute into it 
after the averaging. In this section we demonstrate that the result ( p.25| ) is not satisfactory from 
this point of view. 
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Recall the commutation relations between creation and annihilation parts of the dual fields 

(B) 

[PA(A),g^(/i)] = [Pi,iX),qAilJ')] = \nh{iJ,,X), 

[PdW, qM] = [Pi>Wi Qoif^)] = In h{X, /i), (10.6) 
= ln[/i(A,/i)/i(Ai, A)], 

where h{\, fi) = (A — /x + ic)/ic. It is important that operators pa, Pd, Qa and qn commute with 
each other: 

[p,(A),gfe(^)] =0, j,ke{A,D}, (10.7) 

therefore the vacuum expectation value of any expression containing only the dual fields (pA and 
0D is always trivial: 

(O|^(0^,0Z5)|O) = ^(O,O). (10.8) 

A non-trivial mean value arises if and only if the dual field presents. That is why we paid so 
much attention to this dual field in the previous sections. 

Consider the asymptotic representation ( |10.1| ) from the point of view of its averaging. It is easy 



to see that if we restrict ourselves with the exponential and power law factors, then we immediately 
obtain 

(^O|t^p(<?^A,fc)gt.?"e{0A,<^I3)|O) ^ ^.5^p{0,0)gtJ^,{0,0)^ (10.9) 

since the functionals J-'p and jFg depend only on the fields 0a and 0d- Even if we take into 
consideration the functional JF^^, then nothing will change for the exponential and power law 
factors because of the same reason — JF^^ does not depend on ip. We would like to emphasize 
especially that although the functional J-'c^p is unknown, this fact should not be considered as an 
essential deficiency. The fields 0a and 0d do not infiuence on the exponential and power laws 
factors in the process of averaging. 

At the same time the presence of the functional J-'cii,{4>a, 4>D,'if^ completely changes the situa- 
tion. It depends on the dual field ■0, and it turns out that this operator changes the correlation 
radius, as well as the pre-exponent. Even the simplest factor 6"^^^°^ shifts the arguments of JF^^^, 
J-'p and J-'e due to the evident property 

e^*(^")^(0A(yu), Mf^)) |0) = ^(0A(/i) + Kfx, Ao), 0d(/x) + /i(Ao, /i)) |0). (10.10) 



The contribution of the remaining part of J-'c^ ( 10.5 ) is much more complicated, but it is clear 
that it changes the result (|10.9| ). 



In the present paper we are not going to study the procedure of averaging of the dual fields. This 
is done in 0, where the reader can find the details of the calculation of the vacuum expectation 
value. The main purpose of this paper is to present the asymptotic expression of the operator B in 
the form adjusted to the averaging procedure. The above considerations show that the constant 
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factor J^cip (i-e. the factor, which remains fixed while t goes to infinity) strongly changes the leading 
term of the asymptotics after calculation of the vacuum expectation value. It becomes clear now 
that the leading term of the asymptotics, which had been found in the sections ^ and |^, is not 
sufficient for the description of the correlation function. 

One can ask now, whether we can restrict ourselves with the factor JF^^? Can one guarantee 
that corrections 0{\T?t/t) do not give similar contribution into the vacuum mean value? Let us 
demonstrate that some corrections do provide such non-vanishing contribution. 

In order to find corrections to the expression ( |9.25| ) one can use the differential equations ( p.9| ) 
and (19.101) . Substituting the asymptotic expansions ( |9.11|) , (|9.12|) into these equations we obtain 
an infinite set of the recurrent relations for the operators C,nk and fjnk- The last ones in turn provide 
us with complete asymptotic expansion of the Fredholm determinant and the operator hi2- It is 
clear that corrections, which do not depend on ip, can not change the leading exponential term 
exp{tjFe} of the asymptotics after their averaging. Therefore we need to pay our attention mostly 
to the terms containing the field if). 

Observe that Co(Ao) and r)o(Ao) are proportional to e'^'^'^°^ and e"'^'^'^''^ respectively. This suggests 
the following substitution. 



C(Ao,t) 



C(Ao,t), r)(Ao,t) 



(10.11) 



into (lOI) . We obtain for C, 



dt t 



-0 = 



{u'Y \n\2t) - v" ln(2t) - 2z/V' ln(2t) + ^" + ^'^ ( 



-2(j.'ln(2t)-V^')C' + C"|. 



(10.12) 



Similar equation arises for the operator rj. It is easy to see that the coefficient do of the asymptotic 
expansion ( |9.11|) contains the term 



Cio = ^''(Ao)Co + • • • , 



(10.13) 



and hence 612 contains the term jV''^(Ao)Co- If t goes to infinity, then this term vanishes. However, 
after calculation of the vacuum expectation of ■^/''^/t together with the exponential factor e*'^'= it 
becomes proportional to the positive power of t. Indeed, the action of the operator ip on functions 
of (pA and 0D is equivalent to the differentiating therefore we have 



(0| V'(Ao)e*^'=(^^"^^)|0) ~ te*^^(°'°) + . . . . 



(10.14) 



Thus we see that the term ip''^/t, which was small before averaging and, hence, could be 
considered as a correction, becomes large after calculation of the vacuum expectation. Of course. 
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this term forms only pre-exponent, but does not change the correlation radius. However, it is not 
difficult to see that the higher corrections contain terms ip'^^/t"^, ip'^/t^ etc., which turn into t^, 
etc. after the averaging. Generally, all the terms of type ilj"'/t"^, where n > m, give positive powers 
of t after their averaging together with the exponential factor. The series of positive powers of t 
sums up into the exponent, and as a result the correlation radius changes. 

Thus, we conclude that in order to find correct asymptotic behavior of the correlation function 
it is necessary to have some information about the complete asymptotic expansion of the Fredholm 
determinant and the operator bi2- At least one needs to know the structure of the ip - dependence 
of the higher corrections. The leading term (|9.25 ) itself is not sufficient. This expression correctly 
describes the asymptotics of the operator B, but not the asymptotics of its vacuum expectation 
value. 



11 Modified approach. The 'shifted' saddle point 



We have shown in the previous section that, in fact, in the process of calculation of the vacuum 
expectation the dual field ip is effectively proportional to t. This prompts an idea how to modify 
the method we used for the asymptotic solution of the operator-valued RHP. Recall the jump 
matrix (|3.11|) of the original RHP (a): 



Gn(A) = ?-^(A)Z(A,A)e<^-W|l)(l|; 

Gi2(A) = 27r2(^(A)-l)Z(A,A)e^(^)+^(^)|l)(2|; 



621 (A) 



27r 



7?(A)Z(A,A)e'^^(^)+'^^(^)~^(^)-"(^)|2)(l|; 



G22W = ?-^(A)Z(A,A)e*^W|2)(2|. 

Observe that the field ip enters this matrix only in the combination with the function t(A): ipiX) + 
r(A). All the asymptotic analysis of the operator- valued RHP was based on decomposition of the 
solution in the vicinity of the point Aq = x/2t. This value is the saddle point of the function t(A): 



dr(A) 



0. 



:ii.2) 



A=Ao 



As we have mentioned, one can consider to be effectively proportional to t. Therefore instead 
of point Ao we can consider new saddle point A, which is defined as 



d_ 



-(A) + ^(A)) 



0, 



A=A 



or equivalently 



A = Ao + -V^'(A), 



;il.3) 



:il.4) 
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where prime means derivative with respect to the argument. 

From the point of view of the classical asymptotic analysis this 'shift' of the saddle point does 
not make much sense. In fact, the difference between the original saddle point Aq and shifted one 
A is of order 1/t and it vanishes, when t goes to infinity. The formal solution of the equation (|11.4|) 
can be given by the asymptotic series with respect to 1/t (particular case of ( p.l.6|) ): 



A = Ao + ErTU7 7T^[(^'(Ao)r]. (11.5) 



The replacement of Aq by A means that we have shifted the center of the asymptotic expansion by 
asymptotically small value. Using series (|1 1 . 5|) one can always re-expand any function /(A) into 



the asymptotic series over 1/t with the center at the original point Aq: 

/(A)-/(A„)+|i(^)"^[W'(A„)r/'(Ao)l. (11.6) 



Here we have used the Lagrange's theorem (see, for instance, [Q). As a result, the 'new' asymp- 
totics would of course coincide with the 'old' one. 

Nevertheless this approach appears to be extremely useful when one recalls the role of the 
operator ip in the process of averaging. Indeed, studying the asymptotic solution of the RHP for 
the large t, we automatically take into account the field ip. 

Let us turn back to the RHP (a) and consider the substitution f|4.2| 




X(A) = <I>(A)7^(A), 7^= ^ I . (11.7) 



As before, the operator ^(A) is the solution of the RHP (6), however now we replace Aq by A in 
the corresponding jump matrix: 

g4\) = {e{A - X)Gu + e{X - A){G^2)-') , X e R. (11.8) 

Then one should literally repeat all the constructions of the section ^, replacing everywhere Aq 
by A. Clearly that eventually we shall arrive at the following results for the coefficients of the 
asymptotic expansion of the determinant of g: 

Ao = — / c//isign(A - /i) In (1 - + e*('^) ^is'^('^-^) ) ) , (11.9) 

Ai = / /irf/xsign(A-/i)ln(l-^9(/i)(l + e<^(^)^'sn(A.-A)\\ /^^^^qn 



In comparison with the old results ( |4.38| ), ( |4.39| ) of the section ^ the only difference is that Aq is 
replaced by A. Since A — Aq ~ 1/t, it is clear that new representations for A^- and the old ones are 
asymptotically equivalent. 



60 



It should be noted that the point A must not be necesserely real, so that the jump contour in 
(|11.8|) might be slightly deformed to the complex A with the natural re-definition of the ^-functions. 



Also, the formulae ( |11.9| ) and ( |11.10[ ) should be generally understood as (see ( |1.27| )) 



and 



— r dfiln fl - ^(fi)(l + e'^(^))) , 

TTZ J A ^ ^ ^ ^ 

Ai = ^y'^ /ici/iln(l-?9(/i)(l + e-'^('^))) (11.12) 

-— / /irf/xln(l-^(/x)(l + e^('^) , 

with the integrations go possibly in the complex plane. In all the following formulae we will always 
assume the same interpretation of the integrals involving the sign-functions. 

In order to obtain the Fredholm determinant of the operator / + one need to integrate ( 4.38|) 



and (|4.39|) with respect to x and t respectively. Using 



a. ^^9a, dt^d^-jd^, where = t - ^^"(A), (11.13) 



we have 



~ ~ 1 / \ 

lndet(/ + V^) ^ — / (i/i X - 2/it + iV^'(/i) sign(A - ;u) 
27r j-oo ^ ' 



In (l - ^9(/i)(l + e<^('^)'^'g°(^-^))) . (11.14) 



It is worth mentioning that in the framework of our old approach we obtained the term '?/''(A) 
in the integrand (|11.14|) only after studying the localized RHP and differential equations. As we 



have mentioned already in the previous section, it is this term which is the most important for the 
calculation of the vacuum expectation value, and it is remarkable that the new approach allows 
us to obtain this term already at the first stage of the asymptotic analysis of the RHP. On the 
other hand, one should not be surprised by this fact. The new treatment of the saddle point (|11.3|) 
in fact means that we consider function as it would be proportional to t. Therefor, it is quite 
natural that in addition to the term x — 2/it we have obtained i^l)' . One can also notice that 

z - 2At + #'(A) = i(r(A) + V(A))', (11.15) 



and thus, the integrand in the Eq. (|11.14 ) refiects the fact that the jump matrix in the original 



RHP (a) depends only on the combination r(A) + ^/'(A). Thus, we see that the method proposed 
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in this section permits one to trace the dependency of the Fredholm determinant on the field ifj 
automatically, together with the dependency on the time t. 

We notice again, that in reality the function ip does not depend on t, and it remains fixed when 
t goes to infinity. Hence the representations (|11.14|) and ( [4.41|) are asymptotically equivalent. They 
are not equivalent only when we recall the quantum operator nature of the field ip^X) 



The localized RHP can be considered in the same manner as it had been done previously. We 
would like to mention only one moment. The decomposition of the r(A) + ip{X) at the point A = A 
has the form 

r(A) + ^(A) = r(A) + ^(A) + its{X - Af + 0(A - A)^, (11.16) 

where tg = t — ill)" {A)/ 2 "was introduced in ( |11.13D . We see that variable ts plays the role of the 
new time. Therefore it is quite natural to look for the asymptotic expansion of the solution of 
the RHP with respect to t^, but not to t. Up to replacement t ^ ts and Aq ^ A the solution of 
the new localized RHP formally coincides with the solution of the old one, in particular (cf. ( |9.1|) , 

(11.17) 

(11.18) 

(11.19) 



>12 



7Z/ 



(2t.)— V2e^(^)+^(^)^f2|l)(2|^Jf, 



J2l 



7 



where 7 is equal to 
and 

u = u{A) = 



7 = 2ttZ{A,A){^{A) - l)T{u)e- 



iTTh/ I 3^7^ 
4 



In 



(1 - ^{A)Z{A, A)e'^°(^)) (1 - ^{A)Z{A, A)e^^(^) 



;il.20) 



2m 

/ \ 00 

Of course, now Qq = g^'^\A) and vectors |1), |2) etc. are associated with A, but not with Aq. 

Finally, using the differential equations for 612 and 621 we can find the corrections to the Eqs. 
(|11.17|) , (|11.18|) . The most important is the dependency of these corrections on the function ip. 
Substituting 

b,, = (2t,)-^(^)-i/2e^(^)+^(^)C(A, ts), (11.21) 
b,, = (2t,)'^(^)-i/2e-^(^)-^(^)7)(A,t,), (11.22) 

31), we obtain, for example, for ( 



into equations (B 



-4z-^--(C^C-iz/C) 
dts ts^^'^ 



1 



iy')Hn\2t)-u"\n{2t) C-2i^'C' H'2t)+C' \+0{\n{ts)tf), (11.23) 



where prime now means the derivative with respect to A. Thus, up to the corrections of the 
order we again came back to the equation (|9.9|) . However now we have already taken into account 
the factor e^^^\ and hence we do not need to make a substitution similar to ( 10.11| ). Therefore, 
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the replacement of variables {Xo,t) (A,ts) allows us to remove the dependency on the function 
ip from the terms proportional to tj"^. 

It is easy to see that the asymptotic expansion similar to (|9.11| ) formally solves the Eq. ( |11.23|) 



C = Co(A) + E E L(A)^^^^^, (11.24) 

n=l k=0 '^s 

One can write similar expansion for the operator fj also. 

In order to study the recurrence arising for the coefficients and fj^j^, one needs to know the 
explicit expressions for the terms 0{t~^) in ( |11.23D . It is not difficult to find these terms, however, 
for our goal, it is enough to check that the terms of the order are linear with respect to ip. In 
addition, equation ( |11.23| ) contains the terms of the order tj"^, which are quadratic with respect to 
ip. It is clear that for given n, the coefficients C„fc and r^^^ are polynomials of ip and its derivatives 
of the {n — 1) degree or less. In other words, in the framework of the modified approach any 
large time correction, as it is expected, has the form where m < n. Therefore, as we have 

explained in the previous section, these corrections remain small even after their averaging with 
respect to the auxiliary vacuum, and hence we can drop them out. 

Summarizing the above considerations, we obtain, in the case of negative chemical potential, 
the following asymptotics for the operator B, 

i3 = C(0^,0^,A)(2g-^^e^(^)+-(^) 

xexpj^y" (^x - 2At + z^/''(A)^ sign(A - A) 

X ln{l - ?9(A) (l + e<^W«'g°(^-^)') } ^aI (l + ^(In^ t/t)) . (11.25) 



In distinction to the result (|9.25| ), this asymptotic equation takes into account all the corrections 



which make a contribution into the leading term after the averaging. Corrections to the expression 
11.25|) remain small after calculation of their vacuum expectation value. 



Comparing the results ( |9.25| ) and ( |11.25| ) one can observe that in the last case the constant 



factor C(0D, (pA, A) depends on the shifted saddle point A, which in turn implicitly depends on ip. 
As before this constant factor is not explicitly defined since it contains the operator ^q'^''. In the 

(c) 

framework of the old method Qq depended only on the dual fields (pA and (po, hence this operator 
did not infiuence the leading term of the asymptotics. However now ^q*^^ implicitly depends on 
the field ip, and it might change the asymptotics. Since we do not know explicit expression for 
the operator ^q'^^ we might worry about calculation of its vacuum expectation value. However, 
it is easy to see that the contribution of the constant factor C{(pr),(pA,A) remains constant after 
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averaging. Indeed, remembering the decomposition we have 



C{<p, A) = Ao) + E ^ [Yt) [(^'(Ao))"C"(0, Ao)] . (11.26) 

All the terms of this series behave as {ip/t)"'. Their vacuum expectation values together with the 
exponential factor are some constants: 

(0| e*^'=(<^^'<^«)|0) ~e*^'=(°'°)(C(l) + •••), t-^oo. (11.27) 

Thus the replacement Aq by A in the functional C does not change the exponential and power 
laws of the asymptotics of the vacuum mean value, but changes only the common constant factor. 
Since in this paper we do not analyse this factor, we can simply replace the new saddle point by 
the old one in C. Similar arguments allow us to replace ts by t in the power law factor. However 
we should keep the shifted saddle point A in the exponential factor (i.e. functions sign(A — A) and 
r(A) + iIj{A)), and the power z/(A) since all these objects enter the result together with t and Int 
respectively. 

Our final answer in the case h < can be represented as follows, 
S = C(0^,0^,Ao)(2t)-^^e^(^)+-(^) 

xexpj^y (^x - 2At + z^/^'(A)^ sign(A - A) 

X ln{l - ^9(A) (^1 + e<^W^'s°(^-^)^ } dxj (l + ^(In' t/t)). (11.28) 

In the case of positive chemical potential the situation is quite similar. The asymptotics of B 
is given by the equation 

B = C{<j)D, <I)A, Ao)(2t)-4^e'^(^^)+-(^i) 

X exp 1^ y" (^x- 2Xt + iip'{X)j sign(A - A) 

X ln{l - t9(A) (l + e'^(^) ^'^s^ix-A)^ } dAj (l + 0{1/Vt)) . (11.29) 

The general structure of this result is the same as for negative chemical potential; namely, we have 
the exponential law, the power law, and the constant factors. However, all these factors slightly 
differ from the ones in ( |11.2^ ). First, the integral in the exponent is taken along the contour F (see 



Fig.|]). Secondly, we have the sum r(Ai) + ^/'(Ai) instead of r(A) +'?/'(A). Thirdly, the pre-exponent 
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powers of t also are different. The reasons of these difference were explained in the sections ^ and 
i 

In conclusion we would like to touch briefly the free fermionic limit. In this limit the coupling 
constant c goes to infinity, and one can put all dual fields equal to zero. In this case Eqs. ( |11.28|) , 



11.29) exactly reproduce the results of |29 



Summary 



We study the operator- valued RHP (^.1|), ( |3.11| ). This RHP allows to find the large time and 



long distance asymptotics of the Fredholm determinant, which in turn describes the temperature 
correlation function of the local fields of the QNLS model out off free-fermionic point. Although 
we use the advanced technique based on the nonlinear steepest descent method the main 



idea of the scheme applied in the present paper mostly coincides with the approach considered in 
29| for the free fermionic limit of the model. However, as we saw, the operator-valued RHP is 



essentially more complicated, and not only from the technical point of view. In particular, the 
solution of the RHP was given up to the function q{X,u,v), which is defined as a solution of a 
certain integral equation.. At the same time, in spite of this incompleteness, we were able to obtain 
explicit expressions for the exponential and power laws of the operator B. 

A special problem arises because of the dual fields. Due to the presence of these quantum 
operators one needs to take care not only of the asymptotics of the operator B itself, but also of 
the asymptotics of its vacuum expectation value. This leads us to the non-standard asymptotic 
analysis, which includes the shift of the saddle point. 

The remaining stage in the calculation of the correlation function is the averaging of the 
results ( [I1.28| ), ( |11.29| ) with respect to the auxiliary vacuum. The method for this averaging was 



developed in |^ , although in that paper the contribution of some corrections were not considered. 
Nevertheless, one can definitely state that the calculation of the vacuum expectation value leads to 
nonlinear integral equations, closely related to the equations of the Thermodynamic Bethe Anzats. 
The final result for the asymptotics of the correlation function can be formulated in terms of the 
solutions of these equations. We shall consider these questions in the forthcoming publications. 
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A The Fredholm determinant 



The complete determinant representation for the correlation function ( p. . 1|) has the form 



(vl/(0,0)vl/t(a;,t))^ = -^(0| / " k,{u,v)dudv\0). (A.l) 

27r det (/ - IKt) 



The integral operator I + V acts on the real axis as 

(/ + V) o = /(A) + r V{X, (A.2) 



where /(A) is some trial function. The kernel V{X,fi) is equal to 

1 

A — /i 



V(X,fi) = - / du{E+{X\u)E^{fi\u) - E4X\u)E+{fi\u)). (A.3) 



Here 



^ , 1 Z(m,A) / e--^^^") e-<^°(") , ,,,,, 

E_ (A|n) = —Z{u, X)e-2(M^)+M^)-^W-W) (A.5) 
27r * 

The representation ( [A. 1| ) contains also one more Fredholm determinant, however the last one 



does not depend on time t and distance x, as well as on dual fields: 

2c 

(A -71)2 + c' 



Kt{X, /i) = , ^J mm- (A.6) 



Starting from the kernel V{X,fi) one can construct the jump matrix of the operator-valued 
RHP, using the standard procedure M 



( -E^{X\u)E^{X\v) E^{X\n)E^{X\v) \ 



(A.7) 



V - E_{X\u)E_{X\v) E_{X\u)E+{X\v) ) 
In order to obtain the jump matrix ( 3 .21) , one need to make the following transformation 

G(A)=x°+(A)G™(A)(x°(A))"\ (A.8) 

where 

After this transformation we arrive at the RHP (a), considered in the section]^. 
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B The regularized integral operator 



The regularization of the RHP considered in the section ^ corresponds to new Fredholm deter- 
minant det(J + V^). The kernel K(A,/i) of new integral operator has just the same form as the 
original one 



Here 



K(A,/x) 



1 



A — yU 
1 



du{EUX\u)E'_{^\v) - E'_{\\u)EU^\v)). 



£1(A|..) = -Z(.,A)Z(A,A) 



(B.l) 



[B.2) 



E%{\\u) = E'{\\u)E'^{\\u) 



and 



E'{X\u) 



Z{u,X) 



-MO 



+ 



^- X-iO ^- X + iO 



( -E%{X\u)El{X\v) E%{X\u)E%{X\v) \ 



G*"(A) = 1 + 2711 



\ - Et{X\u)Et{X\v) Et{X\u)EX{X\v) ) 



(B.3) 



(B.4) 



It is easy to see that the limit e — > is well defined in all formulas. In particular — > V . 
The original jump matrix G*" has just the same structure as the matrix ( A. 7 ) 



(B.5) 



In order to obtain the jump matrix one need to make transformation similar to (|A.8| ) with 

matrix 



(B.6) 



The expressions of the logarithmic derivatives of the regularized Fredholm determinant in terms 
of the coefficients h and c ( |2.19| ), ( |2.2(J| ) can be obtained in the same manner, as it had been done 
in The differential equations ( p.21|) also are obvious corollaries of the regularized RHP. 
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